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Abstract. A low-level nocturnal wind maximum is shown to exist over extensive and nearly undisturbed 
rainforest near the central Amazon city of Manaus. Analysis of meteorological data collected during 
the 198S and 1987 Amazon Boundary Layer Experiments (ABLE 2A and 2B) indicates the presence 
of this nocturnal wind maximum during both the wet and dry seasons of the Central Amazon Basin. 
Davtime wind speeds which are characteristically 3-7 ms between 300 and 1000m increase to 
IS m s ' shortlv after sunset. The wind speed maximum is reached in the early evening, with i wind 
speeds remaining high until several hours after sunrise. The nocturnal wind maximum * closely linked 
to a strong low-level inversion formed by radiational cooling of the rainforest canopy. The mght-timc 
inversion extends up to 300 m with strong vertical shear of the horizontal wind below the inversion 
top and uniformly strong horizontal winds above the inversion top. Frictional decoupling of the a. 
above the inversion from the rough forest below, however, is responsible for only part of the observed 
increase Surface and low-level pressure gradients between the undisturbed forest and the large Amazon 
river system and the citv of Manaus are shown to be responsible tor much of the nocturnal wind 
increase. The pressure gradients are interpreted as a function of the thermal differences between 
undisturbed forest and the river/city. The importance of both the frictional decoupling and the horizon- 
tal pressure gradient suggest that the nocturnal wind maximum does not occur uniformly over all 
Amazonia. We suspect that stronger low-level winds are pervasive under clear skies and strong surface 
cooling and that, in many places <i.e„ near rivers), local pressure gradients enhance the low-level 
nocturnal winds. 


1. Introduction 

The lower atmosphere of the global tropics is generally embedded within the 
equatorial trough and is mainly thought to be a region of low wind speed; note 
the widely-used term ‘doldrums-. It came, therefore, as a considerable surprise to 
these investigators to find low-level, nocturnal wind speeds in the near-equatorial 
central Amazonian rainforests to be in excess of 15 m s 1 during both the wet and 
dry seasons. 

Equatorial continental regions of the earth's surface beyond the immediate 
environs of a few cities and settlements remain relatively inaccessible. The rainfor- 
est. with its attendant size (canopies at 40 to 50 m) and biological and physical 
conditions, presents a hostile environment, particularly in the wet season, for 
scientific measurements. Despite these obstacles, considerable work has been done 
in equatorial rainforests. Most has focussed upon the biological aspects of this 
large and diverse ecosystem. The work that has dealt with the atmosphere has 
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been mainly directed at the larger scale dynamics. In the Amazon Basin, for 
example, several researchers (Kousky, 1979: Kousky and Ferreira, 1981; Kousky 
and Kagano, 1981; Paegle. 1987; Salati and Vose, 1984: Silva Dias et al., 1987; 
Nobre et al 1987) have been responsible for pioneering work. The British- 
Brazilian experiments of the early 1980's (Shuttleworth et al., 1984, 1985) were 
the first in Amazonia to address the detail of the hydrologic cycle and the microme- 
teorology of the rainforest. 

In 1985 and 1987, the U.S. National Aeronautics and Space Administration 
(NASA), as part of its Global Tropospheric Experiment (GTE), joined with the 
Brazilian Instituto de Pesquisas Espaciais (INPE) and the Instituto Nacional de 
Pesquisas da Amazonia (INPA) to mount two major experiments in central Ama- 
zonia. These experiments, called the Amazon Boundary Layer Experiments 
(ABLE), were conducted in the dry (July- August 1985) and wet (April-May 
1987) seasons and were directed at establishing the base atmospheric chemistry 
state over the forest and wetlands of the Amazon basin (Harriss et al., 1988, 1990; 
Garstang et al., 1990). We shall draw upon these experiments to document the 
structure and nature of the nocturnal velocity field in the boundary layer over the 
rainforest of central Amazonia. 

The occurrence of a pronounced nocturnal wind maximum over the central 
Amazon basin (CAB) was first noted by Greco et al. ( 1989). Garstang et al. ( 1990) 
and Oliveira ( 1990) have also discussed this feature. In this paper, we shall describe 
the nocturnal evolution ot the boundary-layer winds together with the associated 
temperature and pressure fields. We shall then present a mechanism which ac- 
counts for the night-time wind speed maximum in the first kilometer above the 
forest. 

Nocturnal low-level jets and wind maxima have been documented in many 
locations around the globe (Bonner, 1968; Wipperman, 1973). The low-level jet 
of the Great Plains region of the United States has been extensively studied by 
among others. Hoecker (1963, 1965); Blackadar (1957); Wexler (1961) and Bonner 
(1968). Other mid-to-higher latitude nocturnal low-level jets have been detected 
in southeastern Australia (Clarke et al., 1971), southern England (Thorpe and 
Guymer, 1977). northwestern Germany (Kraus et al., 1985). the United States 
Gulf coast (Hsu, 1979), and the coastal regions of northeastern United States and 
the Canadian Maritimes (Dickerson and Neumann. 1982). In the tropical and 
subtropical regions ot the world, how'ever, only a few documented cases of low- 
level nocturnal jets or wind maxima are in existence. DeSouza et al. ( 1971 ) showed 
remarkable accelerations below 500m over the island of Barbados (near 13° N. 
60° W). These accelerations were not, however, strictly nocturnal. In the Daly 
Waters region of Australia (16°16' S, 133°23' E). Clarke and Brook (1979). Brook 
(1985), and Garratt (1982, 1985) all have documented a nocturnal jet based upon 
measurements from the 1974 Korrin expedition. Virji (1981) found a nocturnal jet 
in the western Amazon basin near the foothills of the Andes. 

Although these strong nocturnal winds form in numerous and varied locations. 
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they have a number of common features. The nocturnal wind maxima all form 
shortly before or after sunset, within the nocturnal planetary boundary layer (often 
at heights between 300-700 m above ground level), and can reach speeds which 
are 2-3 times larger than daytime values. Each low-level jet or wind maximum, 
however, may owe its origin to a different mechanism or set of mechanisms. 

One of the mechanisms most often responsible for nocturnal maxima is the 
inertial oscillation induced when frictional constraint imposed by daytime turbu- 
lence is released by the formation of an inversion near sunset (Blackadar. 1957). 
This allows the wind above the inversion to reach supergeostrophic levels. Other 
prime mechanisms include: winds associated with surface radiational cooling over 
sloping terrain (Lettau, 1967, 1983; Zeman, 1979); advective accelerations (Paegle 
and Rasch, 1973); variations in the synoptic-scale geostrophic wind (Mahrt, 1981); 
and the generation of mesoscale pressure gradients in response to heating and 
cooling over sloping or non-uniform terrain (McNider and Pielke, 1981). 

The ABLE, located in the CAB essentially on the equator, precludes inertial 
oscillation (period too long), Coriolis force (approaches zero for horizontal mo- 
tion) or terrain (large-scale slope ^1:2000 km) as forces driving the nocturnal 
maxima of the low-level wind. However, horizontal surface temperature gradients 
produced by the confluence of the large Negro and Solimoes rivers to form the 
Amazon may play a role in generating the nocturnal wind maxima of the CAB. 
In addition, the effective removal of surface frictional drag by the formation of a 
low-level nocturnal radiation inversion may also be responsible for the increase in 
horizontal wind speed between day and night. 

2. Experimental Design and Data Collection 

A manned tethersonde and radiosonde station was operated in ABLE 2A from 
July 15 to August 7, 1985 at Ducke Forest Reserve, approximately 20 km northeast 
of Manaus, Brazil (Figures 1 and 2). Vertical profiles of temperature, horizontal 
wind speed and direction, pressure, and humidity up to 1000 m were taken by a 
sonde-equipped tethered balloon. The ascent rate of the balloon was 1ms and 
measurements of horizontal velocity were made every 2 m and those of tempera- 
ture, humidity and pressure every 4 m. Each sounding lasted between 15 and 
20 min and both up and down profiles were taken. Three to four vertical soundings 
were made during most hours. In addition to the tethersondes, rawinsondes using a 
high-resolution Atmospheric Instrumentation Research (AIR) sonde and a digital 
recording and processing system were launched every three hours. These rawin- 
sondes ascended to 500-400 mb and provided measurements of wind speed and 
direction, temperature, and specific humidity. Hours of operation were mainly 
confined to the daytime, except for six days when measurements were taken 
throughout the 24 h. The data from these six days will be a main focus in this 
paper. 

In the ABLE 2B wet season experiment, from April 1 to May 14, 1987, Ducke 
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Fig. 1. Location of the ABLE 2A and 2B experiments within the central Amazon basin as indicated 

by the darkened triangle. 



Fig. 2. Orientation, dimensions, and location relative to the Amazon river system of the mesoscale 
triangle network (blackened triangle in Figure 1) employed during ABLE 2B. 


Reserve was one of three sonde-equipped stations forming a mesoscale triangle 
of approximately 1000 km 2 (Figure 2). The mesoscale network consisted of 3 
stations (Ducke, Embrapa, and Carapaha) located at the triangle corners and one 
station (ZF-1) located in the center. The tethered balloons ascended up to 500 m 
at Carapaha and Embrapa and up to 1000 m at Ducke. No balloons were launched 
from ZF-1. Vertical profiles of the boundary layer during ABLE 2B were not 
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obtained as often as during the earlier dry season experiment. At most, 2-3 
vertical profiles were taken during each hour of observation and there were many 
periods of up to 2 hours without a vertical tethersonde profile. Reasons for this 
included logistics, different experimental objectives, and preventive maintenance 
on strained equipment. Rawinsonde launches occurred 4 times a day at Ducke 
and Carapana (0800, 1100, 1400, and 1700 LST) and 6 times a day at Embrapa 
(2000, 0200. 0800, 1100, 1400, and 1700 LST). 

The four stations shown in Figure 2 were also equipped with Portable Auto- 
mated Mesonet (PAM) systems mounted on 45 m towers, placing them 5 m above 
the forest canopy. Measurements of horizontal wind speed and direction, tempera- 
ture, humidity, and pressure were taken at a frequency of 1 Hz and compiled into 
1 min averages. 

Similar to ABLE 2A, most of the ABLE 2B experiment was conducted during 
the daytime. However, a four-day period of continuous 24-h operation was carried 
out near the end of the experiment to examine the structure and evolution of the 
transition (dusk, dawn) and night-time boundary layer during the wet season. 

3. Observed Boundary-Layer Structure 


3.1 ABLE 2A 

Meteorological conditions during the 1985 ABLE 2A experiment were typical of 
the CAB dry season when undisturbed conditions with suppressed convection tend 
to dominate. Night-time observations in the dry season were taken on 6 occasions: 
July 25-26, July 26-27, July 30-31, July 31- August 1, August 1-2, and August 
2-3. With the exception of the initial night-time observation period (July 25-26), 
the tethersonde data coverage was fairly continuous throughout the night. 

Composite soundings for the ABLE 2A experiment were created from the 3- 
hourly rawinsonde launches made at Ducke Reserve. Figure 3 displays the compos- 
ite wind speed profiles below 700 mb for soundings taken at 2000, 0200, 0800 and 
1400 LST.* Between midafternoon (1400 LST) and early evening (2000 LST), 
there is a noticeable increase in wind speed throughout the lower troposphere. 
The strongest increase is found below 1000 m/900 mb where winds increase from 
4_5 to 7-8 m s“ 1 with a nose-like feature found between 950-900 mb at 2000 LST. 
After this initial acceleration, the winds below 1000 m remain relatively constant 
through most of the night before decelerating after sunrise. 

The composite wind profiles from the rawinsonde shown in Figure 3 offer some 
evidence of a boundary-layer nocturnal low-level jet or wind maximum (WM 
hereafter). However, the rawinsonde speeds are averaged over a depth of 200- 
300 m and give a smoothed profile of the boundary-layer winds, which perhaps 

* Since sunrise and sunset in the ABLE mesoscale triangle region (2-3 °S latitude) is near 6 a.m. and 
6 p.m. during the entire year, and the subject of this paper is closely tied to the diurnal cycle, we shall 
use Local Standard Time (60th meridian) throughout. 
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WIND SPEED (M/S) 

Fig. 3. Composite rawinsonde wind speed profiles measured during ABLE 2A at Ducke Reserve for 
0800 LST (dash-dot), 1400 LST (dotted), 2000 LST (solid) and 0200 LST (dashed). 

does not capture the full extent of the nocturnal wind speeds. Table I shows a 
comparison between rawinsonde and tethersonde measured wind speeds at similar 
heights and near-coincident times. During the afternoon (1400 LST), there is little 
difference between the values of wind speed measured by tethersondes and by 
rawinsondes. This is mainly an indication of the well-mixed afternoon boundary 
layer. At night (0200 LST), however, the rawinsondes, unlike the finer resolution 
tethersonde, are unable to capture the strong nocturnal winds and their vertical 
shear. This is illustrated by the 2-5 m s -i difference between the tethersonde and 
rawinsonde wind measurements. The absence of detailed measurements with high 
vertical resolution in the lower atmosphere of the CAB before the ABLE may be 
a reason why a nocturnal WM has not been previously detected. 

However, as seen in the 0200 LST raw insonde composite wind speeds for all 6 
nights and the 3 nights with a strong nocturnal WM (Figure 4), the rawinsondes, 
despite their poor vertical resolution, are still able to capture the nocturnal WM 
on nights when the tethersonde data indicated the existence of a strong maxima. 
On nights when the WM is more modest, the raw'insonde smooths the vertical 
profile to the point where any evidence of a nocturnal WM is lost. To document 
the existence of this WM, we shall concentrate on the more detailed measurements 
taken by the tethered balloons. 
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TABLE I 


Average wind speed values in the 1000-900 mb layer as measured by near- 
coincident tethersondes (7) and rawinsondes ( R ) at Pucke during ABLE 
2A. Several times are given and a mean difference ( T - R) is calculated 



No. of 

rawinsondes 

launched 

Rawinsonde 
(W in ms l ) 

Tethersonde 
(T in m s 1 ) 

T- R 
(ms 1 ) 

14(H) LST 

975 mb 

12 

2.9 

3.9 

1.0 

950 mb 

13 

4.9 

4.4 

-0.5 

925 mb 

13 

4.7 

4.7 

0.1 

900 mb 

11 

5.6 

6.4 

0.8 

2000 LST 

975 mb 

4 

3.5 

8.0 

4.6 

950 mb 

6 

7.0 

8.6 

1.5 

925 mb 

6 

7.8 

9.1 

1.2 

900 mb 

3 

3.7 

8.8 

5.6 

0200 LST 

975 mb 

5 

2.4 

8.2 

5.8 

950 mb 

6 

7.3 

9.5 

2.2 

925 mb 

6 

8.0 

10.5 

2.5 

900 mb 

2 

6.6 

8.5 

1.9 

0800 LST 

975 mb 

8 

1.8 

5.6 

3.9 

950 mb 

11 

4.1 

6.6 

2.6 

925 mb 

10 

5.4 

7.9 

2.6 

900 mb 

9 

5.0 

9.1 

3.6 


Tables Ha— b summarize the salient features of the diurnal evolution of the 
boundary-layer wind and temperature tethersonde profiles for the 5 cases (omitting 
the August 2-3 synoptic event) of day-to-night operations. The nocturnal WM 
and a surface-based temperature inversion exist on all 5 nights, although with 
varying intensities. In addition, there is in all 5 cases a common wind direction 
shift from SE/ESE to NE/ENE between day and night. The wind shift is especially 
strong during the 3 nights with substantial increase of wind speed from day to 
night (July 30-31, July 31-August 1, August 1-2). 

The night-time feature illustrated in Table II and Figures 3 and 4 is examined 
in greater depth by concentrating on 2 cases of a strong nocturnal WM. 

After a one-day suspension of operations, observations resumed at 1800 LST 
on July 30. The formation of a strong surface-based temperature inversion 
(+2.5 °C/150 m) shortly after sunset can be seen in Figure 5a. This inversion lasts 
throughout the night, weakening only slightly by 0600 LST. The top of the inver- 
sion is situated near 380 m above mean sea level (MSL) or 300 m above ground 
level (AGL) (Ducke is 78 m above MSL). This inversion is much stronger than 
the one noted on the partly cloudy nights of July 25-26 and 26-27 due to the 
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WIND SPEED (M/SI 

Fig. 4. 0200 LST composite rawinsonde wind speed profiles during ABLE 2A at Ducke for: all six 

nights (dashed); and the three nights with strong nocturnal wind accelerations (solid). 

more intense radiative cooling produced by the clear skies which were noted by 
weather observations. 

A large increase in horizontal wind speed at and above the inversion top 
accompanies the stable conditions in the lower 300 m. As illustrated bv Figure 5b, 
winds above the top of the inversion increase from 6-12 m s” 1 near 1800 LST 
to 12-15 m s through the night. These winds are uniformly strong up to at least 
1000 m. Below the inversion top, there is a layer of strong vertical wind shear 
with winds increasing from near calm at the surface to 1 1-13 m s 1 at 300 m AGL. 
The nocturnal inversion and WM, along with the separation of layers noted on 
July 30-31 are all stronger than during the nights of July 25-26 and 26-27 when 
broken cloudiness lasted throughout 0200 LST (see Table II). 

Similar to the July 25-26 and 26—27 cases, however, the stronger horizontal 
winds continue through 0900-1000 LST of the next morning (July 31) until the 
inversion breaks down and adiabatic lapse rates resume (Figures 6a, b). After 
1000 LST, turbulent mixing and surface heating once again produce a well mixed 
layer with wind speeds of <7 m s _l . Shortly before sunset, the surface-based 
inversion has once again formed. With the inversion formation, a significant 
increase in the wind speed is seen before 1800 LST with winds reaching 9ms * 
between 250 and 400 m AGL. 
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TABLE II 


Tethersonde-measured boundary-layer temperature (°C) and wind conditions (ms 1 and 
degrees) at Ducke during 5 nights of ABLE 2A operation. Conditions are noted for the 
surface, inversion top (300 m) and 700 m 


7/25-7/26 

7/26-7/27 

7/30-7/31 

7/31-8/1 

8/ 1-8/2 


(a) Temperature 


1600 LST 


sfc 

28.7 

27.5 

25.4 

29.8 

- 

300 m 

27.0 

26.4 

27.3 

28.5 

- 

700 m 

23.4 

22.8 

24.5 

24.7 

- 

0000 LST 

sfc 

22.1 

22.7 

23.9 

22.7 

22.2 

300 m 

23.2 

23.3 

26.4 

26.7 

26.0 

700 m 

23.1 

21.7 

23.7 

23.9 

24.2 

08(H) LST 

sfc 

23.6 

22.5 

23.7 

22.2 

21.3 

300 m 

23.8 

24.0 

26.0 

26.2 

25.0 

700 m 

21.9 

22.3 

23.1 

23.3 

23.8 

1200 LST 

sfc 

26.4 

28.3 

29.4 

- 

30.6 

3(H) m 

24.6 

26.7 

27.7 

- 

28.8 

70(1 m 

21.8 

23.1 

23.9 

- 

24.7 

(b) Wind 
1600 LST 

sfc 

2.4/84 

Ll/106 

1.1/95 

4.2/94 

- 

300 m 

3.8/106 

3.1/70 

7.3/98 

6.3/114 

- 

700 m 

4.7/111 

2.6/25 

12.0/87 

7.0/95 

- 

(HH)O LST 

sfc 

0.4/80 

1.2/50 

2.0/55 

0.5/44 

1.2/10 

300 m 

4.5/- 

4.7/66 

11.7/59 

10.1/90 

7.9/11 

700 m 

9.5/94 

5.8/103 

10.7/92 

10.8/85 

11.7/10 

0800 LST 

sfc 

2.3/45 

2.1/52 

3.4/50 

0,6/80 

0.6/344 

300 m 

6.3/50 

8.4/44 

11.9/74 

10.6/90 

7.4/345 

700 m 

8.5/93 

8.7/63 

12.0/80 

8.8/96 

9.5/0 

1200 LST 

sfc 

4.1/26 

2.6/130 

5.3/120 

- 

4.5/150 

300 m 

5.1/30 

3.8/148 

3.8/48 

- 

7.5/136 

700 m 

7.5/33 

4.1/90 

4.1/90 

- 

7.9/135 


The night of July 31 -August 1 was almost identical to the previous night (Figures 
7a, b). Clear skies and radiation cooling enhance the surface-based temperature 
inversion which extends up to 300-320 m AGL and lasts the entire night. This 
inversion exhibits a lapse rate of +1-27100 m and also weakens slightly by morn- 
ing. Strong nocturnal winds develop above the inversion top with speeds of 10- 
14ms' 1 between 300 and 1000m. Once again, strong stratification is seen with 
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the inversion separating the lower atmosphere into two layers. Strong wind shear 
exists between these two layers. 

3.2. ABLE 2B 

The rainfall of the April-May 1987 wet season was organized on scales ranging 
from the convective to the synoptic, with most of the rain (80%) produced by 
meso- to synoptic-scale systems originating in the eastern Amazon basin and the 
northern Atlantic coast of Brazil (Greco et al ., 1990). However, Greco el al. 
(1990) also show that 40% of the days were dominated by fairly undisturbed 
conditions with only scattered airmass-type convection. These days are similar to 
the ABLE 2A dry season period. With this in mind, a goal of this work was to 
see if the nocturnal WM seen during ABLE 2A also existed during the wet season. 
In addition, with the expansion of the data gathering network to 3 manned stations 
during ABLE 2B, the spatial extent of this phenomenon could be examined. 

The 24 h of continuous observations began at 0600 LST, May 5 and ended at 
1500 LST, May 9. Tethersonde and rawinsonde measurements were taken at all 
three corner stations of the mesoscale network (Figure 2). As an example, we 
shall present a series of tethersonde profiles for May 5-6, a day and night charac- 
terized by generally undisturbed conditions similar to the dry season. 

Figures 8a, b and 9a, b show tethersonde profiles of temperature and wind speed 
at Ducke during the afternoon of May 5 and the night of May 5-6. In the late 
afternoon (1706 LST), the wind speeds through the entire surface-1000 m layer 
art - 3-7 m s . No inversion is present and the near-surface (approximately 20 m) 
temperature is 28 °C. By 2219 LST (Figure 9a), a strong surface-based inversion 
( + 1-2 C/100 m) is in place, extending up to 250 m AGL; wind speeds above the 
inversion top have increased to 8-12 ms -1 (Figure 9b). Although slightly weaker 
than during ABLE 2A, night-time wind speeds during the wet season are twice 
as high as daytime values. Below the inversion top, strong wind shear is also 
present. The inversion and strong low-level winds continue through the night (0432 
LST) before diminishing after sunrise (0712, 0807 LST). These characteristics are 
very similar to those observed during the dry season and are substantiated by 
further wet season cases not shown. 

Examination of the tethersoundings at the other stations in Figure 2 shows that 
the low-level nocturnal WM develops over the entire area. Figures 10a, b illustrate 
the evolution of the surface-based inversion and significant day-to-night wind speed 
increase which occurred at Carapaiia during the evening of May 5. Above the 
inversion top, wind speeds increase to near 11ms *. Due to equipment failure, 
no tethersonde data exist after 0100 LST on May 6. Nocturnal WM similarly' evolve 
over Embrapa but with evidence that the nocturnal maxima at both Carapana and 
Embrapa are not as large as at Ducke. 

The greatest difference seen between the wet and dry season development of 
the nocturnal W M is related to the increase of disturbed weather during the w r et 
season. More frequent and larger amounts of cloudiness reduce nocturnal cooling 
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Fig. 10. Vertical profiles of (a) temperature and (b) wind speed measured by tethersonde at Carapana 
during May 5-6, 1087. Height is in meters above ground level. 


and the intensity of the inversion and thus the associated maxima in the wet 
season. 

Another factor influencing the magnitude of the nocturnal WM is the prevailing 
wind direction. Figures lla-c present a 24-h time series contour of rawinsonde- 
measured wind direction on, respectively, May 5-6, 6-7, and 7-8 at Ducke. The 
nights of May 5-6 and 6-7 were characterized by a strong low-level nocturnal 
WM. On May 5-6, there was a shift from a southerly component before 1700 LST 
to a more northerly component during the night and earlv morning. A shift in 
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Fig. I la. 

Fig. 11. Time-height contours of rawinsonde measured wind direction at Ducke during (a) May 5- 
6. L987, (b) May 6-7, 1987, and (c) May 7-8, 1987. 


boundary-layer winds from ESE to ENE also occurs on May 6-7, as does a 
significant nocturnal WM. In contrast, Figure lie shows that, on a night with a 
relatively weak WM (May 7-8), the predominant wind direction in the boundary 
layer after sunset was 85-120°. Wind directions measured at Carapana and Em- 
brapa reveal the same wind shifts as found at Ducke. 


4. Mechanism for the Observed Nocturnal Wind Maximum 

A necessary condition shown in Section 3 for a strong nocturnal WM is the 
establishment of a surface-based inversion near or after sunset. This inversion, 
which extends up to 300 m AGL, is produced by significant radiative cooling at 
the surface during nights with clear or mostly clear skies. 

With the formation of the inversion, a two-layered fluid is established during 
the night-time hours. Below the inversion top is a layer of strong wind shear and 
mechanical turbulence as winds increase from near calm at the suiface to 10— 
12ms 1 at the inversion top. The layer above the inversion top is effectively 
decoupled from turbulence and surface frictional affects within the inversion layer. 
Such decoupling leads to a relative maxima of the winds above the inversion. 
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The wind speed increase associated with the frictional decoupling is not, however, 
sufficient to account for the observed nocturnal WM. 

The strength of the nocturnal maxima of the wind is also controlled by a 
thermally-driven pressure gradient and associated wind generated between the 
forest sites and the nearby Negro-Solimoes-Amazon river system in the presence 
of the frictional decoupling. The river-land circulation mechanism was also pro- 
posed by Oliviera (1990). Figure 2 shows that Ducke is approximately 20 km north 
of the river system while Carapana is more than 50 km from the Amazon river. 
The rivers reach widths of >15 km near their confluence and the temperature 
difference between the air over the rivers and that over the undisturbed forest is 
potentially large enough to generate an observed mesoscale circulation. 

With only indirect measurements (aircraft) made over the rivers, Ducke-Cara- 
pana is used as a measure of the river-to-forest gradients in surface temperature 
and pressure. Ducke, which lies 20 km north and east of the rivers and north of 
the City of Manaus, will be under the influence of both the rivers and the city. 
Carapana, however, lying 65 km upwind of the city and 45 km north of the Amazon 
is regarded as a representative rainforest site. 

Airborne radiometric temperatures taken during the ABLE 2A and 2B experi- 
ments show the rivers to be 2-3 °C cooler than the forest canopy during the day 
and 6-8 °C warmer at night. Figure 12 shows a 24-h time series of the ABLE 2B 
composite horizontal temperature and pressure differences between Ducke and 
Carapana (D-C) at the surface as measured by PAM towers. During the night 
(1800 LST-0600 LST). it is warmer at Ducke than at Carapana. After 0800 LST, 
the gradient starts decreasing but it is not until 1000-1100 LST when it is warmer 
over the Carapana forest site (C) than at the Ducke river/city site (D). The forest 
station reaches its maximum temperature at 1400 LST and is warmer than Ducke 
until 1700 LST. The similarity in the trends and slight lag (~1 h) between A l and 
IP indicates that the pressure gradient is controlled by horizontal temperature 
differences between the land/forest and river/city at night. 

The pressure gradient (D-C) becomes negative and reaches a maximum 
(-0.25 mb/50 km) during the late afternoon and early evening (1900 LST). This 
maximum in surface horizontal pressure gradient coincides with the increase in 
wind speeds. The negative pressure gradient between Ducke and Carapana implies 
that the wind generated by the pressure gradient force should be directed from 
Carapana (higher pressure) to Ducke (lower pressure). This is indeed the case. 
As illustrated in Figure 1 1. on the two nights with the strongest WM, a wind shift 
occured from ESE to ENE during the period of the intensifying pressure gradient. 
With a wind direction from the ENE, the air will be directed from Carapana to 
Ducke. Figure 12 indicates that, for the remainder of the night and early-to-late 
morning, the pressure gradient, like the temperature gradient, is fairly constant 
^ — 0.08 to -0.14 mb over 50 km). During this period, the nocturnal w-inds shown 
in Section 3 remain strong but with no large increases. Late morning to early 
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Fig. 1_. ABLE 2B daily time-series composite ot the temperature (solid) and pressure (dashed) 
differences between Ducke and Carapana (D-C) as measured by the PAM towers. Pressure and 
temperature values were recorded every minute. 


afternoon sees a reversal in both temperature and pressure gradients along with 
a diminishing of the boundary-layer winds. — 

The horizontal acceleration of the wind can be expressed in the usual fashion: 

dv - 1 _ 

— = — VP-/KXV + F. m “ 

dr p 

For the purpose of this paper, we shall represent the frictional drag within the — , 

tropical PBL in the simplified form 


where "a" is a frictional coefficient with some dependence upon wind speed, 
surface roughness and static stability. Neglecting the Coriolis effect as very small 
near the equator, the equation of motion reduces to: 
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TIME (LSI) 

Fig. 13. Daily time-series of the pressure difference between Ducke and Carapana (D-C) as measured 

by the PAM towers on May 5, 1M87. 


— = — VP-flV. (3) 

d / p 

After a large initial acceleration near sunset, the wind speed throughout the 
PBL was shown above to remain relatively constant during the night. With the 
acceleration near zero at night, Equation (3) reduces to 

— YP= -aX . (4) 

P 

From the PAM-measured surface pressure at Ducke (D) and Carapana (C), it 
was seen that the D-C pressure difference during the late afternoon/early evening 
period of wind increase varied between 0.15 and 0.60 mb, significantly higher than 
seen in the composite (Figure 12). As an example, the PAM-derived pressure 
difference between Ducke and Carapana on May 5 is shown in Figure 13. After 
the early evening WM, A P values remain fairly constant through the next morning 
for all 3 cases, with mean values of 0.15, 0.25 and 0.20 mb on, respectively. May 
5-6, May 6-7 and May 7-8. 

The frictional coefficient "a" can thus be calculated for these cases using Equa- 
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TABLE HI 


Observed pressure differences (D-C) and wind accelerations along with predicted values of total wind 
acceleration, acceleration by the pressure gradient force, and deceleration by friction. Computed using 
the simplified model proposed in Section 4 



May 5-6 

May 6-7 

May 7-8 

Ave. nocturnal 

0.15 mb 

0.25 mb 

0.20 mb 

A vc. nocturnal wind speed 

Urns 1 

1 I 111 s 1 

6ms 1 

a 

2.8 x |() 5 S 1 

4.7 x 10 % 1 

6.9 x 10 % 1 

Observed Day/Night dc/d/ 

5ms 1 in 4 h 

5ms 1 in 4 h 

3ms 1 in 4 h 

Ave. \r during wind speed increase 

0.30 mb 

0.30 mb 

0.20 mb 

V during wind speed increase 

7ms 1 

7.5 tn s 1 

5ms 1 

Accel, by PGF (in 4h) 

9ms 1 

9ms 1 

6ms 1 

Damping by Fric. (in 4h) 

4.7ms ' 1 

5.1ms 1 

3.4 m s ! 

Calculated Day/Night dv/dr 

4.3 ms ‘ in 4 h 

3.9 ms * in 4 h 

2.6 ms in 4 h 


tion (4) and the mean nocturnal values of AP and observed wind speeds. These 
results are summarized in Table III for Ducke. With the mean value of the 
frictional coefficient “ a " from the three nights, a predicted day-to-night increase 
can be calculated using Equation (3) and the observed AP and wind speed during 
the nocturnal maxima period. These values are also summarized in Table III. 

Predicted values of wind speed increases during the period of acceleration on 
all 3 nights are within 1.5ms” 1 of the observed wand increases (see Table III). 
The results suggest that the initial forcing of the WM is due to the increase 
produced by a combination of the thermal/pressure gradients between Carapana 
(forest) and Ducke (river/city), and the decrease in the frictional forcing. Through 
most of the night, little wind increase is observed as the pressure gradient and 
frictional forces are in near-equilibrium. 

After sunrise, a decrease in wind speed responds to the imbalance between the 
frictional and pressure gradient terms. Figures 12 and 13 show' that the pressure 
remains higher at Carapana than at Ducke through midday, well after the observed 
decrease of the wind. This decrease in wind speed between the surface and 1000 m 
is caused by the increase in friction resulting from the destruction of the surface- 
based inversion by strong surface heating and recoupling through turbulence to 
the rough forest. This decrease can be explained by Equation (3) in terms of a 
higher daytime value of the frictional coefficient k 7/" than used during the night. 

As an example, an observed decrease of 2.5 ms 1 was noted at the inversion 
top by Ducke tethersondes between 0627 and 1027 LST on May 7-8. Using a 
value of «=lx 10 4 (slightly higher than the night values) and r = 5.2ms" 1 
and A P = 0.20 mb during the decrease period, the increase by the pressure gradient 
force is 6ms 1 and the frictional decrease is 7.6 ms" 1 . The calculated total 
decrease of 1,6 ms 1 per 4 h is only 1 ms ! below the observed deceleration 
(Table III). 
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5. Summary and Conclusions 

A nocturnal wind maximum in the atmosphere below 1000 m is observed over the 
extensive and nearly undisturbed rainforest north and east of the central Amazon 
city of Manaus. Wind speeds which are characteristically less than 5ms 1 at 3- 
500m above the canopy during the day, increase to 10-15 ms 1 after sunset. A 
maximum is reached during the early evening, after which wind speeds remain 
high until after dawn. 

The nocturnal wind maxima occur during both the dry and wet seasons of central 
Amazonia, and are closely linked to the formation of a strong low-level inversion. 
The inversion is due to radiational cooling of the rainforest canopy and is strongly 
controlled by cloudiness. Thus, nocturnal inversions and associated wind increases 
are more frequent in the dry than in the wet season. The intensity of the nocturnal 
wind maximum is also influenced by the prevailing wind direction. Wind shifts 
observed during the times of wind speed increases suggest that the forcing is 
related to surface pressure gradients. 

Surface pressure gradients are observed and interpreted as a function of surface 
thermal gradients between undisturbed forest and the large river system and the 
City of Manaus. The surface pressure gradients are shown to be responsible for 
part of the observed increases in wind speed between day and night. Frictional 
decoupling of the air above the inversion from the rough forest below must be 
called upon to explain the total increase. 

The dual nature of the nocturnal wind maximum depending upon both the 
existence of a horizontal pressure gradient and frictional decoupling would suggest 
that the phenomenon is not uniformly occurring over the vast reaches of Ama- 
zonia. We suspect, however, that in many places local pressure gradients are 
enhancing the low level nocturnal winds and that, in general, stronger winds in 
the lower atmosphere are pervasive when strong surface nocturnal cooling exists. 
Thus, the phenomenon has seasonal characteristics which relate to cloudiness. We 
also suspect that in the dry season, in particular, strong low-level nocturnal wands 
over wide regions of Amazonia are important in terms of advection of air above 
the surface-based inversion. The finding that biomass burning products were, for 
example, very widely dispersed over the Amazon basin early in the dry season of 
ABLE 2A may be largely due to this process. 
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SUMMARY 


Research under Agreement NCC1-106 during the interim period 
1 October 1991 to 31 May 1992 has continued to ^ se ^he data 
collected during all three ABLE missions (ABLE I, 1984, ABLE Ila, 
1985; and ABLE lib, 1987) . 

The work reported upon in this interim period fol^s the 
description provided in the previous progress report (September 
1991), P i.e., global interactions between the rain forest of the 
Amazon Basin and local and regional processes interior to the 
Basin itself. The progress in each of these two areas during 
this reporting period is summarized below and documented m the 
attached appendices. 

1 . Global Interactions 


a. Aerosols 

A major research result documenting the role played by 
Saharan dust in the Amazon rain forest is now published in the 
journal TELLUS (Swap et al., 1992). The paper is attached in 
Appendix I. The results of this work are summarized bel 

Saharan dust is shown to enter the Central Amazon Basin 
(CAB) in bursts which accompany major wet season ram systems. 

Low level horizontal convergence feeding these rain systems dra 
dust from plumes which have crossed the tropical Atlantic under 
the large-scale circulation fields. Mass exchange of air between 
the surface and 4 km over the eastern Amazon basin is calculated 
using rawinsonde data collected during storm events. Mean 
concentrations o£ dust observed by aircraft over the western 
tropical Atlantic are used to calculate the amount of du 
injected into the Basin. Individual storm events inject some 
480,000 tons of dust into the northeastern Amazon Basin. Storm 
and dust climatology suggest that the annual importation of dust 
is in the order of 13 Mtons . In the northeastern basin, this m y 
amount to as much as 190 kg ha’ 1 yr 1 . Deposition of trace 
species such as phosphate, associated with this dust ranges from 
1-4 kg/ha/yr. Uncertainties in these estimates are not believed 
to be greater than ± 20%. The deposition fluxes f: rom Saharan 
dust are essentially identical to the CAB we ^ d ®P°si ion 
from precipitation in the wet season; a result that implies that 
the major ionic composition of rain water in the CAB wet season 
may be strongly influenced by inputs of material originating 
the African continent nearly 5000 km away. 

The total amount of Saharan dust calculated to enter the 
Amazon Basin is one half to one third of that estimated to cross 
60 ° W longitude between 10° and 25°N latitude. We conclude that 
part of the productivity of the Amazon rain forest is dependent 
upon critical trace elements contained m the soil dust 
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originating in the Sahara/ Sahel . This dependence should be 
reflected by expansions and contractions of the Amazon rain 
forest in direct relationship to expansions and contractions or 
the Sahara/Sahel. Turnover rates for nutrient species deposited 
with Saharan dust in the Amazon Basin suggest a time scale of 500 
to 20,000 years. We believe the dependence of one large 
ecosystem upon another separated by an ocean and coupled by the 
atmosphere to be fundamentally important to any view of how the 
global system functions. Any strategy designed to preserve the 
Amazonian rain forest or any part thereof should equally concern 
itself with the interrelationship between the rain forest, global 
climate and arid zones well removed from Amazonia. 


Research has begun during this reporting period to determine 
the extent to which nutrients contained in aerosols being emitted 
by continental north Africa contribute to the ocean ecosystem. 


b. Energy 

The role of the Amazon Basin as a global center of action 
has been quantified in terms of deep convective transports and 
the contribution of these transports to the planet's hea 
balance. Results of this work are contained in a paper now 
submitted to the MONTHLY WEATHER REVIEW (Appendix II) entitled 
"Amazon Coastal Squall Lines - Part II: Heat and Moisture 

Transports" . 


Mesoscale heat and moisture budgets are presented for three 
synoptic scale Amazon coastal squall lines (ACSL) over the 
central Amazon Basin during the 1987 wet season. The budgets are 
obtained from volumetric analyses of mesoscale rawinsonde data 
from the second Amazon Boundary Layer Experiment (ABLE II ) . 
Classification of cloud and rainfall components are based on data 
from Portable Automated Mesonet (PAM) towers, GOES images, and 
radar. Satellite estimates of the total active cloud area, 
coupled with the mesoscale budget calculations, are used to 
determine the upper level vertical heat transport for an entire 
ACSL as a function of life cycle. 


The instantaneous upper level heat transport in a mature 
ACSL equals nearly 20% of the theoretical export requirement for 
the equatorial trough zone. A little more than half of e 
system-wide transport occurs in the dynamically active anvil 
cloud. Amazon coastal squall lines are shown to exist over the 
basin simultaneously in various stages of formation, maturity and 
decay. When three stages are present, the combined instantaneous 
production of heat exceeds one third of the requirement to 
balance the global heat budget. 

Work has been started on establishing the above findings on 
a scale of the basin itself. To achieve this result we are using 
the large scale cloud fields to identify organized synoptic and 
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mesoscale features. Once identified, then the rainfall and 
diabatic heating associated with these systems will be 
calculated . 

2 . Local and Regional Interactions 


a . 


Nocturnal wind structure 


The paper describing the nocturnal wind fields entitled 
"Low-level Nocturnal Wind Maximum over the Central Amazon Basin 
(Greco et al., 1992, Appendix I) has now been published in the 
journal of BOUNDARY-LAYER METEOROLOGY. 

A low level nocturnal wind maximum is shown to exist over 
extensive and nearly undisturbed rain forest near the cent ral 
Amazon city of Manaus. Analysis of meteorological data collecte 
during the 1985 and 1987 Amazon Boundary Layer Experiments (ABLE 
2a and 2b) indicate the presence of this nocturnal wind maximum 
during both the wet and dry seasons of the Central Amazon Basin 
(CAB) . Daytime wind speeds which are characteristically 
between 300 m and 1000 m increase to 10-15 m s' shortly after 
sunset. The wind speed maximum is reached in the early evening, 
with wind speeds remaining high until several hours after 
sunrise. The nocturnal wind maximum is closely linked to a 
strong low-level inversion formed by radiational cooling of the 
rain forest canopy. The nighttime inversion extends up to 300 
with strong vertical shear of the horizontal wind below the 
inversion top and uniformly strong horizontal winds above the 
inversion top. Frictional decoupling of the air above the 
inversion from the rough forest below, however, is responsible 
for only part of the observed increase. Surface and low-level 
pressure gradients between the undisturbed forest and the large 
Amazon river system and the city of Manaus are shown to be 
responsible for much of the nocturnal wind increase . ' The 

pressure gradients are interpreted as a function of the therm 
differences between undisturbed forest and the river/city. Th 
importance of both the frictional decoupling and the horizontal 
pressure gradient suggest that the nocturnal wind maximum is no 
uniformly occurring over all Amazonia. We suspect that stronger 
low-level winds are pervasive under clear skies and strong 
surface cooling and that, in many places (i.e., near rivers), 
local pressure gradients enhance the low-level nocturnal win 

b. Structure and kinematics of Amazon coastal squall lines 

A manuscript describing the structure and kinematics of 
Amazon coastal squall lines has been submitted to the MONTHLY 
WEATHER REVIEW (Appendix II) . 

This paper describes the meso-to-synoptic scale lines which 
form along the northern coast of South America as seabreeze 
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induced instability lines and which propagate through the Amazon 
basin using data collected during the Apnl-May 1987 ABLE 2b. 

These systems, termed "Amazon Coastal Squall Lines" (ACSL), 
have been noted by others, but details of the structure and 
evolution of the ACSL are limited. The present paper uses 
analyses of GOES satellite, upper air rawinsonde and surface 
Portable Automated Mesonet (PAM) data to describe the structure, 
dynamics and life-cycle evolution of the ACSL. Twelve ACSL were 
sampled during ABLE 2b and three examples are discussed in 
detail . 

The ACSL are discontinuous lines of organized mesoscale 
cloud clusters which propagate across the Central Amazon Basin a 
speeds of 50-60 km h 1 . The ACSL undergo six possible life cycle 
stages: coastal genesis, intensification, maturity, weakening, 

reintensification and dissipation. Analysis also indicates tha 
the ACSL are composed of three distinct components: pre storm 

cumulus (PSC) , leading edge convection (LEC) and multiple 
precipitating cloud layers in the trailing stratiform region 

(TSR) . 

Divergence and vertical velocity calculations indicate deep 
vertical ascent in the LEC and a region of mid-level convergence 
500 mb) inthe TSR. The latter is associated with a weak 
updraft above and an unsaturated downdraft below. The TSR 
vertical motions are an order of magnitude smaller than in the 

LEC. 

Substantial shear in the low-level inflow occurred in all 3 
case studies and, as suggested by model simulations, may play an 
important role in the longevity (24-48 h) of the ACSL. Pre- and 
post-ACSL profiles of 0 C demonstrate that the ACSL stabilize the 
troposphere in their wake and erode a pre-storm minimum 0 e . I 
is hypothesized that the removal of this minimum is accomplished 
not only by direct mixing via vertical motions in the LEC ( hot 
towers"), but also through detrainment in the multiple layered 

TSR. 


c. Gap formation 

Work was begun in this period to document the maximum wind 
speeds at canopy height and relate these winds to the processes 
which cause them and to the results following such winds. The 
creation of gaps or blow-downs in the ram forest is a well known 
phenomena described by ecologists. They attribute a number of 
important processes in the rain forest to gap dynamics including 
diversity and propagation. 
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Saharan dust in the Amazon Basin 
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abstract 

S,n,„„ Pm. a 

.ion tie Ida »« o ... “ Z£Z£L «™. Mean co,ccn.,.,»o.o( dad 

is calculated using rawmsonde data coll ^‘~ , * ^ ^ t0 calculate the amount of dust 

observed by aircraft over the western tropua^ $ome ^qoo tons of dust into the north- 

injected into the Basin, tndividu s suggest that the annual importation of dust 

eastern Amazon Basin. Storm and may amount to as much as 190 kg 

is in the order of U Mtons. In the nor**** abated with this dust ranges from 

and may be as low as ± 20 % The Jj , he w « season; a result that unpl‘«j h ^ 

to the CAB wet deposition flutes f rom CAB wet may be strongly influenced by 

the major ionic composition of rain contin cnl nearly 5000 km away. The total amount 

inputt of mate dal ong.na.mg on e n o 1 12 to 1/3 of ihai etiimated io crosa 

of Saharan dot, =atod.«d » am £ .!.» pad of .he »« 

60" W longitude between 10 and -5 elements contained tn the sod dust 

of the Amazon rain forest by expansions and 

onginaiing in the Sahara/Sa e • 1 relationship to expansions and contractions of 

contractions of the Amazon rain nu . nent specle3 deposited with Saharan dust in the 

the Sahara/Sahel. Turnover rat “ Q / " ^O OOO^rs. We believe the dependence of one large 
Amazon Basin suggest a time scale of500 ^ ltmosph ere to be fundamen- 

ccosystem upon another separated by ^ functions. Any strategy designed to 

tally important to any new of how *< * * ^ equally concern itself with the 

inTeT-reLlmV™w«n a lhe ram forest, global climate and and zones well remov 
Amazonia. 


1. Introduction 

The tropical ram forests of the world are 
generally regarded as large ecosystems m rd»uvely 
delicate nutnent balance. Balance is achieved by 
efficient cycling in a nutrient-limited ^vironment 
(Jordan. 1985; Vitousek and Sanford, 1986). In th 
northeastern and parts of central Amazonia, where 
nutnent poor oxisols and uit.sols are present and 
nutnent supply by nvers is limited l - ° 
airborne soil dust panicles is needed- to achie 
nutnent balance (Jordan, 1985; Salat, and Vose. 


1984; Vitousek and Sanford, 1986). " n ' e am ° un ^ 
of aeolian influx have been bst.mated to be 
relatively high (in the order of millions to ten 
of millions of tons per year). In the case some 
of the critical trace elements th.s mfl ux may be 
crucial to the health of the forest Jordan 1985) 
and Vitousek and Sanford . (1986), for example 
have suggested that phosphorus may be a enttea 
airborne^ element wh.ch could imut annual net 
nnmarY production of the rain forest. 

A substantial body of literature 'dent.fies the 
Sahara and sub-Saharan Sahel as the probable 
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source region of dust found over a large fraction of 
the earth’s surface. Saharan dust transported over 
distances of greater than 5000 km has been shown 
to reach areas in northern Europe (DeAngeiis and 
Gaudichet, 1991; Reiff et al„ 1986) and the middle 
east (Levin et al„ 1980). In addition, there is con- 
siderable evidence mainly provided, by Prospero 
and colleagues that Saharan dust reaches the 
eastern Caribbean (Prospero and Carlson, 1972; 
Talbot etal., 1990), southern North America 
(Savoie and Prospero, 1977) and the northeastern 
coast of South America (Parkin et al„ 1972; 
Prospero et al„ 1981; Prospero et al., 1987; Talbot 
etal., 1990). 

There is also general agreement that transport of 
dust to the far western tropical Atlantic and the 
Caribbean occurs in the northeasterly trades at 
relatively low altitudes, mostly below 6-7 km 
(Prospero and Carlson, 1972; Jaenicke and 
Schiitz, 1978; D’Almeida, 1986; Westphal etal., 
1988; Bergametti etal., 1989). This layer of air. 
called the Saharan Air Layer (SAL), is situated 
most frequently between 850 and 700 mb ( 1.5 and 
3 km) (Prospero and Carlson, 1972; Reiff et aL 
1986; Bergametti et al., 1989). 

Our results suggest that a considerable amount 
of the dust found in the Central Amazon Basin 
(CAB) is of Saharan origin. This Saharan dust is 
transported across the Atlantic under steady state 
(undisturbed) conditions in the 850 to 700 mb 
layer (Prospero and Carlson, 1972). However, 
we suggest that such steady state transport is 
not sufficient to bring the dust into the CAB. 
Synoptic scale weather systems, with associated 
organized low level convergence, are required to 
gather the dust into the Basin. Thus, dust events in 
the CAB are a function of synoptic scale systems 
coinciding with the existence of a dust plume in the 
far western tropical Atlantic. The role played by 
synoptic scale systems in producing pulses of 
Saharan dust in the CAB may also explain why 
Prospero etal. (1981) observed the maximum in 
dust concentrations in Cayenne, French Guiana 
to occur during the height of the wet season. 
Similarly, synoptic scale systems provide a 
mechanism in the far western Atlantic to inject 
dust across an active Intertropical Convergence 
Zone (ITCZ). Penetration of Saharan dust across 
an active ITCZ into the southern hemisphere has 
been deemed unlikely (Junge, 1979) because of wet 
removal of the dust 


Prospero and Carlson (1972) and Prospero 
etal. (1981) have made estimates of the annual 
amount of dust reaching a particular location. 

D* Almeida (1986) has also given estimates of the 
annual amount of particulate matter exported 
northward, westward and southward from the 
Sahara/SahcL Uncertainties, however, remain in 
the estimates of the total amount of dust reaching 
such locations. An associated uncertainty also 
exists in the number and amount of trace elements 
delivered to a given location. 

In the discussion which follows, aircraft 
sampling of aerosols, mass divergences calculated 
from rawinsonde stations and rain and dust storm 
climatology are used to calculate the storm and 
annual amount of dust brought into the Amazon 
Basin. A breakdown of the trace species concentra- 
tions present in the dust is also attempted. The 
resulting annual amounts of dust arc compatible 
with those obtained by Prospero and Carlson 
( 1972) across 60° W between 10 and 25° N latitude. 
The results obtained support the belief that the 
dust observed during the wet season in the central 
Amazon Basin is of Saharan origin. The storm 
convergence mechanism called upon to inject 
dust into the CAB provides a means to calculate 
the event, seasonal and annual amounts of dust 
imported. 

If Saharan dust provides a major source of 
nutrients for the Amazonian rain forests, then a 
series of most important consequences follow. 
Reichholf (1986) has suggested that major fluctua- 
tions in rainfall over northern Africa during and 
following the pleistocene should be reflected by 
parallel expansions and contractions of the 
Amazon rain forest. Petit et al. (1981) have indeed 
reported evidence of increased continental aerosol 
concentrations during the pleistocene, attributed 
to increased size and activity of sand deserts 
located between 30° N and 30° S. Short term varia- 
tions in the amount of particulate matter coming 
off West Africa have been related to the recent 
Sahelian droughts (Prospero and Nees, 1977, 
1986). The possible response of the Amazon rain 
forest to variations in the airborne nutrient supply 
depends upon whether, how, what kind and how 
much dust is imported into the Amazon basin. The 
findings we present below to answer some of these 
questions are based upon a series of National 
Aeronautics and Space Administration (NASA) 
experiments carried out under the Division of 


Anru 

Expc 

2. C 

L 


i 

I 


T 
chcr 
by 7 
imp 
Lay 
Tall 
the 
198 
2b, 
tur’r 
Jut' 
sea 

ant. 

ac: 

50 

sta 

of 

as. 


to 

of 

in 


ca 

as 


di 

a* 

E 

g: 

t: 


Tell us 44B (1992), 2 


SAHAAAN dust in the amazon basin 


135 


^ Pro spero 
of the annual 
ular location. 

dmates of the 
tter exported 
ird from the 
-i\ remain i n 
iust reaching 
ertainty also 
race elements 

■ws, aircraft 
es calculated 
d dust storm 
e storm and 
the Amazon 
s concentra- 
rmpted. The 
compatible 
nd Carlson 
3 N latitude, 
ief that the 
1 the central 
The storm 
n to inject 
to calculate 
nts of dust 

source of 
st s, then a 
ces follow, 
jor fluctua- 
iuring and 
.‘fleeted by 
is of the 
ave indeed 
tal aerosol 
attributed 
id deserts 
erm varia- 
cr coming 
the recent 
:es, 1977, 
azon rain 
rnt supply 
and how 
lasin. The 
sc of these 
National 
(NASA) 
vision of 


Atmospheric Chemistry's Global Tropospheric 
Experiment (GTE). 

2* GTE Atlantic and Amazon Boundary 
Layer Experiment (ABLE) 

The NASA GTE is attempting to document the 
chemical composition of the global troposphere 
byselecting and examining- the atmosphere over 

Layer Experiment (ABLE 1 ) (Gregory et al, 1986 
Talbot et aU 1986) exammed the troposphere o 
the open western tropical aUanac°cean«nJune 
1984. The two Amazon expenments ABLE 2a ana 
2b examined conditions over the reiauve y un is- 
urbed tropical rain forest of the CAB m . Jh. dry 
July-August 1985. and wet, Apnl-May 1987. 
Jons fHamssetaU 1988; 1990L 

Field measurements made in the ABLE 1. - 
and 2b provide information on a range of inter- 
acting sLcs (e.g., space scales from 50 km to 
5000 km) which appear essential to t e un 
standing of both the chemistry and the dynamics 
of this fluid system (Fig. 1). The field «penmen«. 
associated observaaons and results are sum 
martzed in two special issues of the JGR devoted 
To ABLE 2a and 2b (JGR. 1988; 1990). A number 
of these results are used to present a ^oism ° 
introduce Saharan dust into the CAB and to 
calculate the total importation of dust and 

associated trace species. 

Sampling of atmospheric aerosols w« 
ducted from the ground (Artaxo et al.. I > 
at various altitudes (0.2 to 6.0 km) from the NASA 
Eiectra research aircraft (Talbot et al, 1990k The 
ground-based sampling network was located on 
fhree of the towers shown in Fig. l»> at 45 tn, 28 
and 2 m. Sampling duration ranged from ,4 to 
48 h Samples were collected continuously from 
April 8 to May 14. 1987. The concentration of 
particulates and elements contained in soil dust 
(e z St, Al, Fc) was determined on a daily basis. 

Elemental concentrations were determin 
using particle induced X-ray emission (P1XE) 
analysis. Coarse ( > 2 pm) and fine ( < 2 ^mk mode 
particle contributions were determined. Artaxo 
et al. ( 1990) found the following 22 elements m the 
ABLE 2b ground-based aerosol samples: Na, g. 
Al Si P S Cl. X. Ca, Ti. V. Cr. Mn, Fe, Ni. Cu, 
Zn. Br, Rb. Sr, Zr, and Pb. Rain water analysis 
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ra l Basin-scale rawinsonde network (a) and local 
JJe network (b) employed in the Amazon b«,n 
during ABLE 2B. Locauon of the mesoscale ne 
indicated by the darkened mangle in (a). 

was performed on samples collected at 3 of the 
tower locations shown in Fig. lb (Andreae et al, 
1990). Rain water samples were collected on an 

event basis at each station. 

Airborne (NASA Eiectra) sampling of aerosols 
consisted of collection on teflon filters and subse- 
quent analysis of the water-soluble trace species 
using ion chromatography (Talbot « tal ~J 9 V 
The constituents observed included NiUNH,. R. 

Q. NOj, SO,. C 2 0„ CH,SO„ HCOO, and 

CH.COO (Talbot et al, 1990). 

The European Centre for Medium Range 
Weather Forecasting (ECMWF) globa win 
analyses are used for the computation of 10 day 
backward and forward isobanc trajectones. An 
overview of the ECMWF data assm^auon syswn 
may be found in Hollingsworth (1987). The 
ECMWF analyses are available every 6th h with 
horizontal resolution corresponding to 125 km 

19 vertical hybrid levels. 

The 2-D trajectory method is a straightforward 
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Lagrangian following of air parcels in the analyzed 
wind field with wind components interpolated 
linearly in time every 6 h. The wind components 
are further interpolated, both horizontally and 
vertically (pressure), to the actual position of the 
parcel. The interpolated wind components are 
then used to track the parcel forward or backward 
over 6 h in 15-min increments. 

Horizontal wind fields were compiled using 
ECMWF global data together with commer- 
cial aircraft and satellite wind observations 
(T. Knshnamurti, personal communication). The 
data at 850, 500, and 200 mb are used to determine 
the streamline fields over the tropical Atlantic, 
Africa, and South America during the ABLE 2b 
experiment 

Satellite imagery from the European Meteosat 
are used to provide coverage of the tropical 
Atlantic from northern Africa to the coast of 
South America including the arid region of north- 
east Brazil. 

Ground-based rawmsonde measurements, ob- 
tained 4 times daily, were made in a basin scale 
network ( 1000-1500 km apart) shown in Fig. la. 
Local meteorological measurements used in this 
study include rawinsondes from the three comer 
stations of the mesoscale network in Fig. lb 
(approximately 50'km spacing) together with 
rainfall measurements made on the towers at 5 m 
above the canopy. 

Aerosol optical thickness values deduced from 


the Advanced Very High Resolution Radiation 
(AVHRR) satellite measurements (Rao et aL, 
1988) were not available during any of the ABLE 
missions. We refer to post field experiment aerosol 
optical thickness data to confirm the frequency 
and geographic extent of Sahara dust outbreaks. 

3. Analysis and results 

Ground-based measurements of aerosol concen- 
tration and composition during ABLE 2b indicate 
the occurrence of several large pulses containing 
both terrestrial and marine aerosols (Fig. 2a). 
Three maxima were centered around April 12-14 
(experimental days 12-14), April 29-May 1 (exper- 
imental days 29-31) and May 5-8 (experimental 
days 36-40). 

Aerosols associated with these pulses had rela- 
tively high concentrations of alkali and alkaline 
earth elements in the atmospheric mixed layer 
of the CAB (Talbot et al., 1990). Amazonian 
soils, however, are highly weathered and exhibit 
deficits of these elements. This finding suggested a 
non-local source region for the dust, which Talbot 
ct al. (1990) believe to be the Sahara. Pulses of 
Saharan dust were also observed at Barbados 
during ABLE-2b (Talbot et al., 1990). These 
pulses, with concentrations at the surface reaching 
180 /igm' 3 , appeared at Barbados 1 to 2 days 
before Saharan dust was seen in the CAB. 

The temporal variation of major mineral aerosol 
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constituents (e.g^ Fe, AL Si) shown in Fig. 2. are 
compared with the mean regional mesoscale 
(4 station) rainfall in Fig. 2. Both aerosols and 
rainfall occur in distinct pulses. The two major 
pulses in Si coincided with the two major rain 
events. The third weaker aerosol pulse preceded a 
weaker rainfall event by about 4$ h. Other- rain 
events show little or no association with increases 
in dust. On days with little or no rain, however, 
very low concentrations of dust are observed. 

Analysis of the rainfall during ABLE 2b (Greco 
et aU 1990; Garstang eta L, 1990; Swap, 1990) 
showed that more than 80% of the rain was due to 
large (linear scales of 1000 km or larger) organized 
synoptic scale systems. 

Wet season rainwater chemistry associated with 
large systems compared to small or local rain 
systems is different The highest concentrations of 
SOj~ and Cl~ as well as higher pH values occur 
in large organized rain systems compared to local 
rain showers (Andreae et aU 1990). The presence 
of marine constituents (e.g^ Na* and Cl ) in 
these large rain systems suggest that they are 
ingesting low level air from over the western tropi- 
cal Atlantic. Higher pH values in rain water, 
however, have also been found to be associated 
with Sahara dust (Loye-Pilot et aU 1986). There- 
fore, marine chemical constituents and continental 
components from African regions south of the 
Sahara may contribute to the rainwater composi- 
tion oflargc rain systems in the CAB. 

High concentrations of mineral aerosols occur- 
ring in close association with large ram events 
requires an explanation because wet removal 
would suggest the opposite result We propose 
below that low level convergence covering a large 
area and preceding such large rain systems by 
24 to 72 h is responsible for increasing the dust 
concentrations in near coincidence with the 
production of rain from these weather systems. 
This explanation requires the presence of a 
Saharan dust plume in the far western tropical 
Atlantic at the time that the synoptic scale weather 
system is organizing on the northeastern South 
American coast. First, we show from trajectory 
analysis that Saharan dust can reach the eastern 
and central Amazon Basin. Secondly, we calculate 
mass convergence in the low troposphere pre- 
ceding the weather system and will show that this 
convergence can concentrate the dust from the 
preexisting plume residing over the western tropi- 


cal Atlantic Finally, we show that measured 
values of concentrations in the oceanic dust plume 
allow the calculation of both the total mass 
of dust and the mass of the associated chemical 
constituents transported into the Amazon Basin. 

3.L Streamline and trajectory analysis 
The large subtropical anticyclones are indi- 
cators of the nature of the lower tropospheric 
fields of motion linking western Africa and eastern 
South America. The position of the West African 
Subtropical High (WASH) with ridging over the 
central tropical Atlantic has been shown to 
influence low level flow across the Atlantic into the 
Amazon Basin (Talbot et aL, 1990). Fig. 3 shows 
that, with the WASH centered in a prescribed area 
of West Africa, flow at 500 mb extends across the 
Atlantic and into the Amazon Basin. Each of the 
three concentrated pulses of soil dust into the CAB 
shown in Fig. 2 followed a 4-7 day period of the 
WASH so established over West Africa (Talbot 
et aL, 1990). Meteosat visible images (not shown) 
confirm that, during such periods of dominance by 
the WASH, fair weather extends across the low 
latitude north Atlantic. Under a slowly changing 
large-scale anticyclonic field, the streamlines 
approach trajectories. The absence of deep con- 
vection or other forms of strong vertical motion 
along the path provide conditions which are near 
optimal for the calculation of trajectories. 

Two-dimensional horizontal isobaric trajec- 
tories were calculated at 975, 950, 925, 900, 850, 
800, 700, and 600 mb based upon the 6 h ECMWF 
global wind analysis. The 2-D trajectory calcula- 
tions at multiple levels over 10 days were used to 
produce an envelope to bound the possible path 
taken and as an indication of the approximate 
origin or end point of the air parcels. Visible 
satellite imagery and the large scale wind fields in 
the 5—7 day period prior to a soil dust event in the 
CAB suggest near steady state conditions with 
little or no deep convection in the region of dust 
transport (Talbot et aL, 1990). In the absence of 
deep convective clouds, vertical motions in the 
lower troposphere are confined to large scale 
subsidence in the order of 400 m day 1 ( Richl. 
1954; Riehl ct aL, 1951). Forward trajectories were 
calculated from two locations: 

1 5° N latitude 10° E longitude in the Sahel, and 
20° S latitude 15° E longitude in the Namib. 
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LONGITUDE 

fl g i 500 mb streamlines on 17 April 1987 illustrating the West African Subtropical High and associated wind flow 
across the Atlantic Ocean and into the Amazon Basin. 


Backward trajectories were calculated from 

0 latitude 60° W longitude in the CAB. 

Three 10-day periods were chosen, two to coincide 
with the two major dust events seen in the CAB 
and one to coincide with a period of very low dust 
concentrations in the CAB (Fig. 2). The results of 
these calculations are discussed below. 

SJ.L Dust event of 12 April- 1 4 April. Back- 
ward trajectories starting on 12 April 1987 (Figs. 
4a, b) and 14 April 1987 (Figs. 4c, d) and forward 
trajectories starting on 4 April 1987 (Fig. 5) 
suggest that air parcels between 300 and 700 mb 
(Figs. 4b, 4d and 5b) originate from a region near 
Lake Chad. Air below 800 mb is traced back to the 
eastern Atlantic and the West African coast (Figs. 
4a, c). Air below 800 mb originating (Fig. 5a) near 
Lake Chad does not leave the west coast of Africa 
after 10 days of travel time. Air above 700 mb is 
either not traced back to Africa (Fig. 4b) or does 
not reach the Sahel in the backward trajectory 
calculation. In the 10-day forward calculation 
(Fig. 5b), air at 600 mb heads towards, but does 
not reach the mouth of the Amazon. 

Backward and forward trajectories were also 
calculated for the dust event of 29 April- 1 May 
1987 and for a low dust concentration day, 
26 April 1987. We discuss our findings from 
these trajectory calculations without showing the 
graphical results. The backward trajectories 
starting on 29 April 1987 and l May 1987 and 


forward trajectories starting oq 21 Apnl 1987 show 
results similar to the April 12-14 event Backward 
trajectories at 800 to 700 mb once again end up 
near Lake Chad. Above 800 mb and below 700 mb, 
however, parcels show a diverse set of path- 
ways. Forward trajectories with the calculation 
beginning on 21 Apnl 1987 just west and north of 
Lake Chad do not show any connection to the 
Amazon Basin. Backward trajectory results suggest 
that the point of origin chosen for the calculation 
is too far north and east. Trajectories for the low 
dust case show no connection to the Amazon 
Basin. 

10-day forward trajectories were also calculated 
for three separate dates: 4, 16 and 21 Apnl 1987 to 
determine whether dust could be coming from the 
southern hemisphere deserts of the Namib and 
Kalahari. No such connection was found. 

While the trajectory calculations earned out as 
part of this paper do not rule out any other source 
of soil dust, we believe that the results show 
support for the contention that the Sahara/Sahcl 
is the primary source. We find no evidence in the 
trajectory calculations that would support a South 
American location as a source region in the 
Apnl/May time of year. None of the backward 
trajectories calculated from the CAB point 
(0° latitude, 60° W longitude) approach the dry 
zone of northeastern Brazil. Transport from the 
Namib/TCalahari, as well as from NE Brazil, into 
the Amazon basin may well, however, occur 
during astral spring and summer. 
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Fig. 5. 10-day forward trajectories beginning on 4 Apnl 1987. 


3.2. Horizontal mass divergence calculations 

The large triangle in Fig. la shows the relative 
location of three of the six basin scale rawinsonde 
network stations operated during ABLE 2b and 
designed to describe large-scale meteorological 
conditions over the entire Amazon Basin. The line 
from Boa Vista to Belem is nearly parallel to the 
northeast coast of South America. This line is also 
approximately normal to low level inflow into the 
Amazon Basin under conditions which favor both 
the transport of dust across the Atlantic and the 
development of synoptic scale weather systems in 
the Amazon Basin. The Belem, Boa Vista and 
Embrapa triangle over the northeastern Amazon 
Basin is thus positioned to capture any enhanced 
mass inflow into the Basin, drawing lower 
tropospheric air from over the western tropical 
Atlantic. 


Rawinsonde data from the basin-scale stations 
were used to calculate horizontal divergence for 
the Belem, Boa Vista and Embrapa triangle by 
employing a line integral method utilizing the wind 
components normal to each leg of the triangle 
(Brummer. 1978). The horizontal divergence 
calculated for the tnangle serves to determine the 
amount of air and, with measured aerosol concen- 
trations, the amount of dust imported into the 
tnangle volume as a function of the circulation oi 
the synoptic scale storm system. 

Fig. 6 shows the experiment time senes ol 
horizontal divergence at 850 mb and 700 mb. The 
rawinsonde network was not fully functional until 
15 April 1987. The dust event of 29 Apnl-1 May 
1987 is captured by this time senes and is used 
below to calculate total mass inflow into the 
volume above the triangle. The divergence signal 
in Fig. 6 shows a major perturbation about 
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Fir. 6. Time senes of horizontal divergence at 850 mb (a) 
and 700 mb (b) calculated for the Belem-Boa Vista* 
Embrapa triangle using raw insonde (solid) and 
ECMWF (dashed) data. 


the 29 April- 1 May event. Increasing low level 
convergence precedes the event, followed by 
divergence which is coincident with the major rain 
pulses. 

Divergence calculations for the triangle were 
repeated based solely upon ECMWF data. The 
results of this calculation are also shown in Fig. 6. 
While the ECMWF data captures much of the 
temporal variation in the signal, divergence is an 
order of magnitude lower than that obtained from 
the calculations based on the soundings. 

Divergence calculated for the smaller mesoscale 
triangle shown in Fig. lb (Scala et al., 1990) 
yield divergence values an order of magnitude 
higher (10 _ '*s' 1 ) than those obtained from the 


large scale triangle network soundings. Such a 
hierarcbal increase in divergence (10"‘, 10 
10“' 4 s _l ) proceeding from the .large scale flow 
fields, as represented by the ECMWF, to the syn- 
optic fields large scale triangle (Fig. la) to the 
mesoscale (mesoscale triangle soundings. Fig. lb) 
is consistent with other work in the tropical 
atmosphere (see for example the GARP Atlantic 
Tropical Experiment results such, as Frank, 1978). 
We, therefore, use a 48-h mean convergence of 
2.0 x 10 -5 s " 1 as representative of the mas s inflow 
into the volume bounded by the surface to 4 km 
and the sides Boa Vista to Belem to Manaus. 

J J. Calculation of mass transports of Saharan dust 
into the Amazon Basin 

The area of the large triangle in Fig. la is 
approximately 4.0 x 10" m 2 . The volume of air 
from the surface to 4 km over the triangle is 
1.6 x 10" m 1 . If the mean divergence over 48 h 
preceding the dust event is —2.0 x 10 -5 s"‘, then 
the amount of air exchanged during this time is 
( — 2.0 x 10 _1 a — 1 )( 1.728 x 10 1 s) = 3.5 times the 
original volume. In contrast, with large-scale flow 
in the absence of an organized synoptic scale 
weather system (divergence =10 * s 1 ), the 
exchange rate would be an order of magnitude less, 
or only 1/3 of the low level volume would be 
exchanged in 48 h. Under synoptic scale forcing, 
the volume of air contained below 4 km over the 
triangle is exchanged 3.5 times in 48 h equalling 
(3.5 )( 1.6 x 10" m 3 ) = 5.6 x 10" m 3 of air. 

Under northeasterly inflow into the Amazon 
Basin, we may assume that all of this exchange 

occurs across the side of the triangle parallel to 

the coast Fig. 7 shows the dust plume laying off 
l the NE coast of South America. Fig. 8 shows a 
r Saharan dust plume immediately east of Barbados 
i (60°W, 13°N) as observed by downward looking 
airborne lidar (Browell, personal communication). 

: The plume lies between 1.5 and 4.2 km. In the 

: presence of such a dust plume, the air being 

i. brought into the basin-scale triangle by the storm 
c system will contain large amounts of dust. Talbot 
i et al. (1986) showed that these dust plumes 

i contained up to 400 jig m- J , with a mean value 

of 51 jigm -3 for thc surfacc t0 3 - 71cm layer A 
e weighted mean for the layer between 1.5 and 4 km 
i) from Prospero and Carlson ( 1972) gives a value of 

e 45 M g m -3 . We take a mean value of 50 jig m for 

e the dust concentration in such a plume between 
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Fig. 7. Aerosol optical depth (t) on 21 April 1988 showing a dust plume between western Africa and the northeastern 
coast of Brazil. Values for r derived at NOAA/MHSD1S using satellite radiance values obtained from the AVHRR. 
Contour intervals are multiplied by 100. 


the surface and 4 km over a continuous penod of 
48 h. Using this mean concentration, we calculate 
the mass of dust entering the volume in 48 h to be: 

(the volume of air exchanged in 48 h) 

(mean concentration) 

(5.6 x I0 i3 m 3 )(5.0 x 10“ 3 g m _3 ) 

= 2.8 x 10 11 g or 
280,000 tons per storm system. 

A given storm may have a convergence field 
which acts on the dust plume for longer than 48 h. 
However, this may be counteracted by the fact that 
a mean concentration of 50/igm~ J may be high 
for the entire duration of the inflow. Thus, we feel 
that 48 h and 50 /ig m " 3 are representative of the 
average dust event. 

The calculation above assumes the dust enters 
across the side of the triangle (Fig. la) parallel to 
the coast. In fact, the dust will enter across the 
entire NE coast of South America. Aerosol optical 
depth data from the Advanced Very High Resolu- 
tion Radiometer (AVHRR) presented in (Fig. 7) 
shows that the length of coastline across which the 
dust is transported extends from Sao Luis, Brazil 
(44°30'W longitude, 2°30'S latitude) to the Gulf 
of Pana (63°0'W longitude, 10°30'N latitude), a 


total distance of 2500 km. The linear distance of 
the triangle side to the coastline is 1:1.7. Wc thus 
increase our above estimate by 1.7: 

(2.8 x 10 3 )(1.7) 

= 4.8 x 10 3 tons/storm system. 

Greco et al. (1990) have shown that 80% of the 
CAB rainfall during the wet season is produced by 
organized synoptic scale weather systems. The 
average amount of rainfall in the CAB during the 
6 months (December-May) of the CAB wet season 
coincident with Saharan dust outbreaks is about 
1600 mm (Villa Nova et al. f 1976; Marques-Fiiho 
et al., 1981). Swap (1990) shows that the average 
amount of rain per synoptic scale system is 20 mm. 
Thus, approximately 60 synoptic scale weather 
systems occur in the 6-month wet season or ap- 
proximately 10 systems a month. Other evidence 
(Madden and Julian, 1972) show a 5-day periodi- 
city in the equatorial tropics or 6 systems a month. 

The wet season for the eastern Amazon basin 
(EAB) is approximately 2 months longer than the 
CAB wet season. The EAB wet season lasts from 
December to July (Villa Nova et al., 1976), and 
with approximately 2000 mm of rain falling during 
this period. If 80% of the EAB rain is also 
delivered by organized synoptic scale systems, 
approximately 80 storms occur in the 8-month 
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Fig. 8. Lidar- measured aerosol attenuation coefficient near Barbados on 21 July 1984 (courtesy of E. Browell, 
personaJ communication ). The dust plume is seen between 1.5' km and 42 km. 


wet season (10 systems a month). These numbers 
arc supported by Cohen (1989) who reported an 
average annual value of 67.5 synoptic scale 
instability lines in the EAB during these months 
between 1979 and 1986. Cohen’s numbers are 
confined to "instability lines" and do not include 
all synoptic scale systems. 

Aerosol optical thickness measurements made 
from the AVHRR equipped satellite over a 
32-month period from 4 June 1987 to l March 
1990 are used to obtain a measure of the number 
of weeks in which Saharan dust outbreaks are pre- 
sent off the NE coast of South America. Table 1 
shows that 34.1% of the weeks during the 32 
months have aerosol optical thickness vaJues of at 


least 0.3 within the area bounded by 10°N-0°, 
60°W-40°W. Values of 0.3 or greater are repre- 
sentative of heavy aerosol outbreaks (Carlson, 
1979). Table 1 also shows that heavy dust out- 
breaks are seen from December to July. This 
corresponds to the height of the wet season for 
eastern and northeastern Amazon Basin (Villa 
Nova ct aL, 1976) and to the period of increased 
coastal squall line activity (Cohen, 1989). 

Based on the results of rainfall per system and 
the aerosol optical depth analysis, 34.1 % of the 80 
or 27.3 synoptic scale systems have the potential to 
carry dust into the eastern Amazon Basin during 
the extended wet season of December- July. If we 
round off the number of dust bearing storms to 27 
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Tabic l. Vo of weeks with aerosol optical thickness 
values of (a) <0.2, (b)<0.3 and (c)^OJ for the 
area bounded by 10° S and 0° and 60 ° W—40 ° W 
from 4 June 1987 to I March 1990 



a 

b 

c 

January 

27.2 

18.2 

54.6 

February 

0 

16.7 

83J 

March 

9.1 

9.1 

81.8 

April 

0 

14.3 

85.7 

May 

115 

37.5 

50.0 

June 

0 

53.8 

46.2 

July 

7.7 

61.5 

30.8 

August 

66.7 

33J 

0 

September 

83J 

16.7 

0 

October 

83.3 

16.7 

0 

November 

64.3 

35.7 

0 

December 

58.3 

25.0 

16.7 

Overall 

36.2 

29.7 

34.1 


per wet season then the annuaJ amount of Saharan 
dust introduced into the Amazon Basin is 

(4.3 x 10 5 tons per storm) 

x (27 storms/wet season). 

If, all of the Saharan dust introduced into the 
Amazon Basin in a year is accomplished in the wet 
season, then the amount of dust per year is 
13.0 Mtons/year. 

Prospero and Carlson ( 1972) have estimated the 
annual transport of Saharan dust through a wall 
extending from 10° N latitude to 25" N latitude 
along 60° W longitude to range between 25 and 
37 Mtons per year. D’Almeida (1986) has 
estimated annual westward transport immediately 
off West Africa to lie between 45 and 53 Mtons per 
year. Estimates of the total annual West African 
source strength for westward transport vary from 
80 Mtons per year (Prospero and Carlson, 1972) 
to 260 Mtons/year (Jaenicke and Schiitz, 1978) 
with 190 Mtons per year having been reported by 
D’Almeida ( 1986). 

Our results suggest that the amount of dust that 
enters the Amazon basin is between one half and 
one third of the amont of dust found by Prospero 
et al. ( 1981 ) to cross the 60th meridian between 10 
and 25" N. The major uncertainties in our calcula- 
tions are the mean concentration (50 ^gm* ] ) and 
the number of synoptic scale systems (27) which 


bring dust into the Amazon Basin. Mean concen- 
trations of dust in the outbreaks coinciding with 
synoptic scale systems in the Amazon Basin are 
unlikely to drop below -25 j/g m “ J or to rise above 
75 m Extreme variations in mean wet season 
rainfall in the- CAB reach ±20% (Kousky and 
Kagano, 1981). The presence or absence of a single 
synoptic scale systenr in each of the 8 wet season 
months can account for such a variation in rainfalL 
Thus, a maximum or minimum possible variation 
in dust supply due to a variation in storm systems 
is ±8 systems. These extremes, if coincident, give 
a range of 4.4 Mtpns/year to 24.4 Mtons/year. We 
believe that such an extreme range is unrealistic 
For example, all storms added or subtracted will 
not coincide with dust events; nor are concentra- 
tions of dust likely to go down in dry years or 
up in wet years. Both effects reduce the extremes. 
A more realistic range is half the extreme range or 
8.5 to 18.5 Mtons/year. 

3.4. Transport of trace species with Saharan dust 

The mass concentrations of selected water- 
soluble species associated with Saharan dust over 
the western tropical Atlantic and the CAB are 
given in Table 2. Many of the nutnent species 
associated with Saharan dust (e.g., K, NH 4 , N0 3 , 
P0 4 ) are believed critical to the physiology of 
tropical ram forests (Vitousek and Sanford, 1986). 
The mass % contributions of the various species 
range from tenths of a percent to several percent. 

Mass concentrations of Na and Cl are enhanced 
over the CAB compared to their levels in the 
Sahara air layer over the western tropical Atlantic 
(Table 2). This apparent enrichment is presumably 
due to entrainment of sea salt into dust-containing 
air masses that are advected inland from the South 
American coast. The sea salt concentration over 
the CAB is on the order of 0.5 /jgm" 3 during 
periods influenced by the injection of marine air 
(Talbot et ai., 1990). Besides Na and Cl, NH 4 and 
S0 4 appear to be enriched 2-fold and NO, 3-fold 
over the CAB compared to the Sahara air layer 
over the western tropical Atlantic We attribute 
these differences to inputs of nitrogen-rich biomass 
burning emissions along coastal Africa (Del mas, 
1982). The levels of NH 4 , N0 3 , and SO* 
associated with Saharan dust over the tropical 
Atlantic appear to be strongly influenced by 
African biomass burning products during the 
winter-to-spring tune of the year (Savoie et al.. 
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Table 1 Mass concentrations (% ) of selected trace 
species associated with Saharan dust over the 
Western Tropical Atlantic and the Central Amazon 


Basin 

Specie* Average (%) 

(a) Western tropical Atlantic *' 


Na 

K 

NH 4 

ci 

no 3 

so 4 

po* 


0.37 

0.17 

0.18 

0.32 

0.51 

1.29 

0.037 


(b) Central Amazon basin 

Na 

l.l 

K. 

0.29 

nh 4 

0.44 

Cl 

1.9 

NO, 

1.6 

so 4 

2.4 

po 4 

N/A 


Range (%) 


0.24 -0.85 
0.10 -0J3 
0.01 -0.27 
0.19 -0.50 
0.21 -1-5 * 
0J5 -3.0 
0.017-0.071 


0.42-1.8 
0.12-0.46 
0.08-0.91 
1.3 -2.6 
0.20-16 
0.22-3.7 
N/A 


•» From Talbot et a L (1986). 
** From Talbot ct aL ( 1990). 


of other reported Saharan dust depositional rates; 
13kgha~ l yr~ l at Miami* Florida (Prospero 
etaL, 1987) and 77kgha-'yr-‘ in the French 
Alps (DeAngelis and Gaudichet, 199 1 ). 

Using the mass concentrations of selected trace 
species associated with Saharan dust over the CAB 
(Table 2) and our above estimate of dust deposi- 
tion to the Amazonian ecosystem, we find that 
most species are deposited at a rate of 1—4 kg 
h a -> yr _l (Table 3). These deposition fluxes are 
essentially identical to the average wet deposition 
fluxes measured at two stations within the triangle 
area during the ABLE 2B experiment (Table 3). 
This result suggests that precipitation scavenging 
is the principal removal mechanism of Saharan 
dust associated material from the ■ atmosphere 
during the’ CAB wet season. During the CAB wet 
season, the average residence time of aerosols in 
boundary layer air relative to precipitation 
washout is 0.25-1 day (Andreae et aL, 1990). Our 
analysis suggests that the major ionic composition 
of precipitation (e.g, NH„ N0 3 , SO,) tn the CAB 
wet season may be strongly influenced by inputs of 
material originating on the African continent 
nearly 5000 km away. This transcontinental bio- 


1989). Our ABLE l measurements over the western 

tropical Atlantic were conducted during summer- 
time and exhibit minimal influence from African 
biomass combustion products; a dramatic sea- 
sonal signal in this source appears in the annual 
record of aerosol composition at Barbados (Savoie 
et aL, 1989). Thus, Saharan air intrusions into the 
CAB during the wet season appear to reflect 
combined chemical inputs from mineral aerosol, 
sea salt, and possibly biomass burning products. 

The annual rates of supply of trace species were 
estimated assuming that all the dust entering the 
Boa Vista/Manaus/Belem triangle is deposited 
within that area. In assessing this assumption, it 
should be noted that Manaus lies more than 
1300 km inland from the NE coast of South 
America. If ail the Saharan dust introduced into 
the triangle (Fig lb) in a year was effectively 
deposited in that area, its deposition flux would be: 

13 x 10 9 kg yr — _ (90 ha' ' yr -1 . 

(l.7)(4.0x 10 7 ha) 

This value, 190 kg ha ' 1 yr ' \ is on the high end 


Table 3: Comparison of deposition fluxes 

(kgha~‘ yr~‘) of selected trace species associated 
with intrusions of Saharan dust and precipitation in 
the Amazonian wet season (December- May) 


Dust Intrusions*' Precipitation*’ 

Species average range average range 


Na 

11 

0.80 -3.4 

1.7 

0.73-29 

K 

0.55 

0.23 -0.87 

1.3 

0.41-23 

nh 4 

0.84 

0.15 -1.7 

0.71 

0.38-1.4 

Cl 

3.6 

25 —4.9 

19 

1.4 -5.1 

NO, 

3.0 

0.38 -4.9 

1.4 

0.90-3.7 

so 4 

4.6 

0.42 -7.0 

1.8 

1.2 -3.1 

po 4 

0.07 

0.032-0.14*' 

0.04* * 

0.01-0.12' 


*’ Calculated using mass concentrations reported in 

m 1 Based on measurements reported in Table 1 from 
ndreae ct aL ( 1990). 

«> Based on upper and lower qua rule concentrations 
werved at two sites within the triangle area (Andreae 
aL, 1990k 

«> Es tima ted using avenge mass concentrauo 
able 2 ) reported by Talbot etaL (1986k 
•’Based on unpublished data from AHLt -o 
Talbot and M. Andreae, personal communication. 
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geochemical link appears to be severed during 
the dry season when the large-scale meteoro logical 
regime which favors transport of Saharan air 
masses into the CAB migrates to the north. The 
aerosol and precipitation chemistries in the CAB 
dry season support such a scenario (Talbot et al-, 
1988; Andreae etaU 1988). Long-range transport 
of aerosol-associated species to the CAB appears 
to explain the higher wet season precipitation 
fluxes of Na, Cl, and K. compared to the dry 
season. It also seems to be responsible for the 
similar seasonal precipitation fluxes of NH*, NO,,, 
and SO* despite a five-fold contrast in rainfall 
amounts. Our estimate of long-range chemical 
inputs to the CAB are significantly lower than 
ReichholTs (1986), who suggests at least 2 orders 
of magnitude higher dust transport to the western 
tropical Atlantic atmosphere than any reported 
estimate. 

Table 4 compares our deposition fluxes to 
reported annual internal fluxes (the amount of 
nutrients internally recycled) and standing 
stock values (the total ecosystem amount of a 
nutrient) for the Amazon Bastn. Inputs of K 
and P associated with Saharan dust aerosols 
represent from a few tenths to 6 % of the annual 
nutrient turnover due to litterfall, and from a few 
thousandths to a few tenths of a percent of the total 
ecosystem standing stock. The nutrient influx of 
tiese and other trace elements into the EAB and 
CAB does not appear important on short time 
scales (10‘-10 2 years). However, with respect to 
the long-term geochemical mass balance of the 
Basin. Saharan dust and associated aerosols are 
important as an outside source of elements that 
help replenish those lost from the system. This 
importance is amplified in a tightly constrained 
nutrient system such as the EAB and CAB. 


The time required for the imported K and P to 
be equal to the amount of their respective 
Amazonian standing stocks is estimated using 
our deposition fluxes to evaluate the potential 
biogeochemical importance of Saharan dust inputs 
to the Amazon rainforest Dividing the standing 
stock of K in the Amazon oxisol rainforest 
-380 kg ha -1 (Jordan, 1985), by our maximum 
and minimum values of K influx (Table 3), we 
arrive at values ranging from 440 to 1700 years. 
The same calculation for P, using a standing stock 
value of 220 kg ha" 1 for the CAB (Goiley et aU 
1975) results in values ranging from 4700 years to 
20,000 years. These calculations do not account for 
losses from the system, such as the Amazon River 
outflow (Gibbs, 1972). Hydrologic loss values 
(Franken and Lcopoldo, 1984) of P from the CAB 
vary from an order of magnitude less than our 
deposition values to the same order of magnitude. 

These results suggest that if Saharan dust 
aerosol is a principal source of these trace species 
to the EAB and CAB, then the effect on the 
biogeochemical cycling would be on the long term 
(10M0 3 years). Other possible mechanisms may 
play additional roles in this seemingly storm event 
driven system. The accumulation of nutrients in 
rain water stored in epthytes may be one such 
mechanism. Increased tree falls in tropical rain 
forests are noted during the wet season (Oldcman. 
1972; Brokaw, 1985). Urge trees heavily loaded 
with epiphytes may therefore be prune candidates 
to fall during the CAB wet season. Tree fall is 
most likely to occur in the synoptic scale storms 
called upon to gather the dust into the Basin. 
Downdrafts from deep convective clouds in these 
storm systems reach speeds of 20 to 25 m s 
(45-55 mph). The creation of gaps in the rain 
forest at the time of maximum soil dust influx may 




Deposition*! 

input 

(kg ha -1 yf 1 ) 

0.011-0.047 
0.23 -0.87 

41 Proctor ( 1984). 

KJinge (1976). 
el Goiley et al. (1975). 


Internal flux 
from litterfall 
(kg ha" 1 yr " 1 ) 

1.4— 4.1*' 
13 -21" 


Input aa */. of 
litterfall 
(min— max) 

0.27-3.4 
U -6.7 


Standing 
stock 
(kg ha "‘) 


Input as 
% of 

standing stock 

0.005-0.021 

0.061-0-23 
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be an important component in the utilization of 
small amounts of trace elements. 

4. Summary and conclusions 

We conclude that Saharan dust penetrates the 
eastern Amazon Basin and is detected in the 
Central Basin. The dust is delivered in intermittent 
pulses which occur in association with major 
pulses in rainfall. Since these large rain events are 
due to organized storm systems with horizontal 
scales of up to a few thousand kilometers and time 
scales of days, we conclude that part of the 
mechanism driving the storm system is responsible 
for the accumulation of dust. 

Low-level (up to 4 km) horizontal negative 
divergence associated with the storm system is 
shown to be sufficient to draw dust into the 
Amazon Basin. Previous work together with 
results provided here show that Saharan dust 
easily reaches the western tropical north Atlantic 
Mean flow into the Basin can neither concentrate 
the dust nor deliver it in the aperiodic fashion 
observed. Storm systems are shown to be suf- 
ficiently frequent to draw upon existing dust 
plumes. The presence of Saharan dust oiT the 
northeast coast of South America is also shown to 
coincide with the high frequency of storms. Not all 
rain events in the CAB are accompanied by dust, 
but all of the dust events we observed occurred in 
association with rain events. 

Observations of the honzontaJ extent, vertical 
depth and concentrations of dust permit an 
estimate of the amount of dust introduced into the 
Amazon Basin in a given storm event. Storm 
climatology yields an estimate of 13 Mtons/year of 
dust introduced into the Amazon Basin. In the 
northeastern basin, this may amount to as much as 
190 kg ha " * yr ~ 1 of dust. The amount of trace 
species supplied with this dust ranges from a few 
kg ha " 1 yr ~ 1 to less than 1 kg ha “ 1 yr “ l . 

The total amount of Saharan dust calculated to 
enter the Amazon Basin is one half to one third of 
that estimated by others to cross 60° W longitude 
between 10° and 25° N latitude (Prospero and 
Carlson, 1972). It is two orders of magnitude less 
than is implied by Reichholf (1986) who has 
estimated quantities of trace species needed by the 
rain forest on nutrient poor soils. This large dis- 
crepancy may be due to an underestimate in the 
total amount of soil dust supplied to the Amazon 


Basin if, as we have, one confines the calculation to 
the Sahara as the only source. The discrepancy 
may also be due to an overestimation of the 
necessary nutrient influxes and a simultaneous 
underestimation of the efficiency in recycling 
nutrients in this ecosystem. 

Whatever the resolution of the above paradox 
is, we conclude that part of the Amazon rain forest 
is dependent upon soil dust which originates in the 
Sahara/Sahel region of western Africa. We concur 
with conclusions drawn by others that such a 
dependence must be reflected by expansions and 
contractions of the Amazon rain forest in direct 
relationship to expansions and contractions of the 
Sahara/Sahel. We hypothesize that such fluctua- 
tions in the size of the rain forest are on dme scales 
which lie between about 500 and 20,000 years. 
Evidence for such a coupled response must be 
present in the records of the pleistocene pluvial 
periods. Such a dependence of one large ecosystem 
upon another separated by an ocean is of funda- 
mental importance to our view of the role of these 
ecosystems in the global atmosphere. A view of the 
role of Amazonian rain forest in the global system 
must take this coupled process into account. Any 
strategy designed to preserve the Amazonian ram 
forest or any pan thereof should equally concern 
itself with the interrelationship between the rain 
forest, global climate and acid zones well removed 
from Amazonia. 
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DESTRUCTION OF THE RAIN FOREST AND 
CLIMATE CHANGE 
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INTRODUCTION 

No a priori consequences can be inferred from the possible destruction of the rain forest 
and the impact of that act upon the global climate. The reason why we cannot project 
such a response is that we have little or no understanding of how the undisturbed forest 
influences the global system. If we do not know how the present system works, how can 
we predict the effect of a changed state of one component on the whole? 

Instead the possible effect of the destruction of the rain forest should be posed as a 
series of questions. Our ability to answer these questions will be a measure of how well 
we may be able to assess the possible global effect of deforestation. 

Three questions will be addressed in this paper. First, are there grounds for assuming 
that the tropical rain forests play a role in the global atmospheric system? Second, are 
the tropical rain forests being destroyed and at what rate? Third, depending upon the 
answers to the first two questions, how do the current rain forests function in the global 
system? A follow on to this last question is then, how well can we predict a response 
of global climate to a change in the rain forest? 

SOME BASIC FACTS ABOUT THE TROPICAL RAIN FORESTS 

While rain forests occur in both temperate and tropical latitudes, the tropical rain forests 
are by far the most extensive. Tropical rain forests occupy about 7% of the land surface 
of the Earth. Nearly half of the world’s tropical rain forests are in Brazil’s Amazon Basin 
covering some 336 Mha. One’s view of the potential effect of an ecosystem occupying 
only 7% of the earth’s surface is greatly changed when it is recognized that that 
ecosystem contains 60% of the total global biomass. The primary productivity of the 
tropical rain forest is some 5 Btons of dry material per year, fixing 6.58 tons of carbon 
per year (Salati, 1987). 

The rain forests are the most diverse of all ecosystems on the planet and are thought 
to contain as many as 30 million species of which only 3 million are known. The rain 
forests contain 2500 tall tree species compared to a few dozen tall tree species in 



temperate zone forests. The abundance of life in the tropical rain forest means that this 
ecosystem produces 1 5-20% of trace gases found in the earth’s atmosphere (Harriss et 
al., 1988). Carbon dioxide, carbon monoxide, hydroxyls, heavy and light hydrocarbons, 
methane and nitrogen species are all produced and are all directly or indirectly involved 
in capture and emission of terrestrial radiation. Primary and secondary aerosols, mostly 
organic in composition, influence the transmission (scattering) of solar radiation and the 
formation of clouds which in turn influence the albedo. Water in these environments is 
in abundance (precipitation exceeding evaporation). High water vapor content in the 
equatorial atmosphere implies high latent heat content. High precipitation (in the order 
of 2 to 3000 mm/year) means that large amounts of latent heat are being converted to 
sensible and ultimately geopotential energy. Fresh water runoff is high, the Amazon river 
alone accounting for 25% of the world’s river discharge. 

THE PHYSICAL SYSTEM 

The equatorial regions which contain the bulk of the tropical rain forests are the firebox 
of the global atmosphere (Riehl and Malkus, 1958). Vertical injection and subsequent 
poleward transport of heat to achieve global thermal balance does not occur uniformly 
around the girdle of the planet. "Hot spots" coincident with the large tropical rain forests 
are found over the Amazon, the Congo (Zaire) and over the warm ocean continent of 
Indonesia-New Guinea. Hadley and Walker cells represent mean circulations which are 
fueled by the intermittent action of deep convection coupled to the forest-atmosphere 
interface at the top of the canopy some 40 m above the ground. The rain forest canopy 
has an albedo near 1 0% allowing 90% of incident solar radiation to be absorbed. Solar 
radiant energy at the surface is used in photosynthesis. Up to 65% of heat may be 
returned to the atmosphere via evapotranspiration; the remainder returned by 
evaporation and sensible heat transfer. Some 20% of the rain falling upon the forest is 
captured, much in unique species such as bromeliads which amount to 10% of the 
biomass of the Amazon rain forest. It is estimated, but not known accurately, that 50 /o 
of the rain falling in the Amazon basin is recycled through the vegetation and surface 
waters. 

The sum total of the physical processes outlined above leads to the conclusion that the 
tropical rain forest is coupled in an intimate way with the overlying atmosphere. The 
overlying atmosphere in turn is connected through the larger scale Hadley-Walker 
circulations to the global atmosphere. 

POLITICAL AND ECONOMIC FORCES 

The rain forests are subject to large political and economic forces. Forces acting on the 
Amazon rain forest illustrate the situation. Much, if not all. of the stress imposed upon 
the rain forest stems from the increase in human population. The population of Brazil 
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in 1990 is near 150 million. By the year 2050 that population is expected to double. In 
1975 it was estimated that 0.56 ha of arable land was required to support each person 
in the Brazilian population. Increases in population and the need to improve nutrition 
(upping the amount of arable land to 0.6 ha/person) means that by the year 2050 an 
additional 180 Mha of arable land will be needed. If such a demand is met by clearing 
the Amazon rain forest it would mean consuming more than half the current area of rain 
forest. 

Brazil is under serious energy pressure. Energy needs are expected to double in the 
next decade. Brazil has no reserves of fossil fuel. Hydroelectric power and wood 
(charcoal) are the most immediate sources of energy in Brazil. The Amazon basin is 
estimated to contain 1000 GW of hydroelectricity. It is expected that 22 GW will be 
tapped by 2000. Dams needed to supply these 22 GW of hydroelectricity would 
inundate 100 Mha of forest. 

Projected use of charcoal, primarily to fuel iron ore smelters in the Amazon basin, is 
estimated to require some 90 Mha of forest over the next 10 years. Unless the efforts 
to successfully employ tree species that can sustain such a demand, this use will rapidly 
encroach upon the existing forest. 

Global biomass burning is currently estimated to produce 4 Bton C/yr" 1 . This compares 
to 5 Btons C yr 1 produced by burning fossil fuel. While not all of the biomass burning 
is due to consuming the tropical rain forests, a large fraction comes from this source. 
In Brazil it is estimated that only 5% of the rain forest timber is extracted and effectively 
used. Ninety-five percent is burnt or left to decompose. 

Estimates of the rate of deforestation range from 0.5 to 3.0% per year with some 
estimates as high as 6% per year. Most studies (Logan et at., 1981; Setzer and 
Pereira, 1986; Salati, 1987; Watson et a/., 1990) suggest that deforestation is 
increasing at an exponential rate. Evidence from satellite imaging suggests that in the 
Brazilian state of Rondonia this may have been the case from 1982 to 1987. There is 
now some possibility that the rate of destruction has slowed. 

The above suggests that the first two questions posed are both answered in the 
affirmative. There is substantial basis to argue that the tropical rain forests are coupled 
to the global atmospheric system and that the tropical rain forests themselves are 
undergoing rapid change. 

The remainder of the paper will address the third question: how do the tropical rain 
forests function in the global system? 
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CHEMICAL AND BIOLOGICAL STATES 


Trace gases like carbon dioxide (C0 2 ), carbon monoxide (CO), methane (CH 4 ), isoprene 
(C 5 H 8 ) and nitrous oxide (N 2 0) are injected into the atmospheric boundary layer 
immediately above the forest canopy under undisturbed and disturbed (burning) 
conditions. 

In the undisturbed state, trace gases generated in the rain forest are well mixed vertically 
by active moist convection. Figure 1 shows the vertical distribution of a trace gas (CO) 
with its source at the surface (forest) under undisturbed conditions and in the wake of 
widespread biomass burning. When burning is taking place, the background 
concentrations of CO are seen to treble and marked layers of trace gas are formed in 
the lower troposphere. 



Fig. 1. Vertical distribution of carbon monoxide in the central Amazon basin under wet (July 12, 1985) and 
dry (August 12, 1985) conditions. Biomass burning is responsible for the enhanced concentrations observed 
on August 12. 


Horizontal distributions of the trace gases under undisturbed and disturbed conditions 
have been shown to be large (JGR, 1990). Measurements in the dry season (burning 
season) show that products of biomass burning extend over the entire Amazon basin 
over linear distances of up to 2500 km. 

The presence of soil dust was also noted in the Central Amazon Basin (CAB) (Talbot et 
at., 1990). Figure 2 shows pulses of soil dust elements (Fe, Al, Si) occurring at or nearly 
at the same time that major rain events are detected in the CAB. 
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Fig. 2. Time-series over approximately 6 weeks in the central Amazon basin of ground-based 
concentrations of terrestrial elements (Fe, A I, SI) together with rainfall. 

PHYSICAL PROCESSES 

The chemical measurements show that gases and particles originating in the forest 
canopy are injected into the atmosphere. The sections below will examine aspects of 
the vertical transport through three atmospheric layers: 

* the canopy-atmosphere interface (100 m deep); 

* the atmospheric boundary layer (2 km deep); and 

* the cloud layer (15 km deep). 

CANOPY-ATMOSPHERE INTERFACE 

Figure 3 shows the diurnal flux across the canopy of sensible heat (w'T'), latent heat 
(w’q’), solar radiation (Rs) and the net radiation (Rn = Rs - outgoing longwave radiation). 
Air temperature (T) and specific humidity (q) are also shown. 

Strong diurnal forcing is noted with latent heat flux reaching 400 W m 2 and sensible heat 
flux 100 W m' 2 soon after midday. Latent heat fluxes decrease to very small values at 
night while sensible heat flux reverses direction with heat being lost by the atmosphere 
to the forest. 

The maximum latent heat flux over the forest is twice the amount exchanged by the 
tropical ocean to the atmosphere. The corresponding sensible heat flux of the forest is 
six times that occurring over the open tropical ocean. The combined latent and sensible 
heat input equals the net radiation and approaches the total solar input. 
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Fig. 3. Time series under fair weather conditions in the central Amazon basin at 5 m above the forest 
canopy of heat flux w'T, temperature (T), moisture flux w’q\ specific humidity (q), incoming solar 
radiation flux Rs, net radiational flux Rn, and rainfall. SR and SS are approximate sunrise and sunset. 


The flux of latent heat during the day is such that water vapor divergence occurs at the 
top of the canopy with specific humidities dropping 2 to 3 g/kg. This "draw-down" of the 
canopy water vapor must be compensated for by nighttime and early morning recharge. 

Other measurements suggest that the greater fraction of the latent energy input is by 
evapotranspiration (65%) as compared to evaporation from water surfaces (35%). The 
large energy injection via water vapor during the central hours of the day is clearly an 
important characteristic of the rain forest. It is likely to trigger a significant buoyant 
response in the lower atmosphere leading to convective cloud development and further 
and deeper vertical transport of quantities out of the canopy. 

ATMOSPHERIC BOUNDARY LAYER 

Figure 4 shows the growth of the mixed layer over the rain forest. For a 5-hour period 
beginning 1 hour after sunrise, the mixed layer grows at a rate of 8-10 cm s' 1 . 
Entrainment across the top of the mixed layer approaches 600 W m 2 which is a 
contributing factor to the drying observed in the surface layer above the canopy in Figure 
3 above. 

Mixed layer growth ceases at about 1 600 LST decaying rapidly to a shallow (< 300 m) 
night-time mixed layer. A fossil mixed layer (well mixed but no longer actively mixing) 
remains. In the absence of cloud convection (particularly precipitating), the fossil mixed 
layer builds concentrations of species on a day-to-day basis. Moist convection acts as 
a siphon drawing down these concentrations. 

The vertical transport of quantities out of the forest canopy depends not only upon the 
fluxes across the top of the canopy but upon the depth and rate of growth of the mixed 
layer. 
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COMPOSITE 



SR TIME (LST) SS 

Fig. 4. A composite of 7 fair days in the CAB showing the mean mixed layer height (solid line) 
and individual daytime mixed layer heights (crosses and squares). The nocturnal mixed layer 
heights are marked with an asterisk. Contours are virtual potential temperature v in °K. 

Mixed layer growth rates as predicted by models (Tennekes, 1973; Betts, 1973; Lilly, 
1 968; Fitzjarrald, 1 982) show a good agreement in the early hours of growth after which 
model produced mixed layer heights are too large. Adjustment to the excessive model 
heights can be achieved by calling upon convective cloud induced subsidence in the 
order of 2 to 3 cm s' 1 . Thus, cloud transports which take over from the mixed layer 
transports are both upward and downward. 

CLOUD LAYER 

Cloud transports are calculated using observations of clouds captured within an 
observational triangle (Figure 5). Rawinsonde measurements were made at the corners 
of a triangle 50 km on a side. Convergence into the triangle was calculated and vertical 
motions in the cloud cells deduced from the fraction of active cloud (10% of total area) 
present. It was assumed that the calculated upward motion in the triangle in the 
presence of clouds occurred within the cloud cores. Vertical motions of up to 5 m s 1 
were obtained in the clouds. 

The observations of cloud transports were then used to help verify cloud model 
calculations of transports. A two-dimensional moist cloud model was used (Garstang et 
al„ 1988; Scala et al., 1990). Figure 6 shows the cloud model results. Figure 6 (a) 
shows a single cell cloud which was found to develop in the very moist environment of 
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DI5TRNCEIKM) 

Fig. 5. Observational triangle 50 km on a side in the central Amazon basin with vertical 
soundings of wind, temperature and humidity on the corners and radar coverage. Radar 
echoes are shown with the cores of the echoes shaded. 

the wet season. Figure 6 (b) shows a multicell cloud which occurred in the dry season 
and intermittently in the wet season. Significantly different vertical transport is achieved 
by the single and multiple cell clouds. Backward trajectories from 12 km elevation in the 
single cell cloud (Figure 7 (a)) show that 73% of the air reaching 12 km originates above 
6 km. Only 15% of the air reaching 12 km originates in the mixed layer. The single cell 
wet clouds do not effectively couple the deeper atmosphere to the surface and mixed 
layers. These clouds, however, are very effective mixers of the middle and upper 
troposphere. 

The multicell cloud occurring in drier conditions (Figure 7 (b)) however, forms a very 
different link between the surface and deeper atmosphere. Ninety-two percent of the air 
reaching 12 km originates below 4 km, 60% below 2 km in the multicell cloud. 

Dogma which has persisted for 30 years since the first formulation of the hot tower 
concept (Riehl and Malkus, 1958) has held that all deep vertical transport in the 
equatorial trough region is performed by the multicell clouds. We believe from our 
findings that the single cell clouds which predominate in wet season conditions create 
a condition that permits mean circulations such as the Hadley and Walker, to perform 
deep vertical transport. 
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Fig. 6a. Model generated single cell cloud showing the boundary of the cloud in terms of 
specific humidity (g/kg). Areas of upward and downward motion are shown by light and 
stippling respectively. Model grid units (abscissa) equal .05 km. 
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Fig. 6b. As for Fig. 6a except for a multicell cloud with potential temperature (oK) instead of specific 
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Fig. 7b. Model generated backward 
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Thus, over the rain forest at least two mechanisms exist for deep vertical transport: 

* multicell clouds occurring in drier unstable atmospheres, and 

* single cell clouds occurring in wetter more stable atmospheres. 

When the latter exists, mean circulations can transport surface species through the depth 
of the troposphere. 
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HORIZONTAL TRANSPORTS 


The physical processes described above have concentrated on vertical transports in the 
atmosphere over the rain forest. Our results suggest that at least one long range 
horizontal transport phenomenon is of importance. Figure 2 showed pulses of soil dust 
coinciding with the major rain events in the Amazon basin. The occurrence of high dust 
concentrations under the rainiest conditions would appear to be a paradox in the sense 
that rain cleanses the atmosphere of aerosols. We explain this paradox by showing that 
the large rain systems of the Amazon basin draw dust into the basin from dust plumes 
that have originated in the Sahara-Sahel and spread across the Atlantic (Swap et ai, 
1991). 

Synoptic scale (~ 1000 km) convergence between the surface and 4 km is needed to 
draw the dust into the basin. Calculations show that some 460,000 tons of Saharan dust 
are injected into the basin by a single synoptic scale storm system. Storm and dust 
climatology suggest that more than 12 Mtons of Saharan dust are introduced to the basin 
each year. In the northeastern basin this may amount to as much as 180 kg/ha/yr. The 
amount of trace elements, such as phosphate, supplied in this dust range from a few 
kg/ha/yr to less than a kg/ha/yr. The delicate energy balance existing in the rain forests 
which have developed on nutrient poor soils suggests that the critical trace elements 
supplied by Sahara dust may be a limiting factor in the survival of the rain forest. 
Dependence of one ecosystem upon an ecosystem on another continent separated by 
an ocean and linked by the atmosphere is a real possibility. 

CONCLUSIONS 

The rain forests inject significant amounts of trace gases into the atmosphere. These 
gases are detected through much of the troposphere and over the region of the rain 
forest. Major regional enhancements of certain gases are noted as a product of biomass 
burning. 

Local, regional and large scale circulations suggest that the rain forests are coupled to 
the global atmosphere. The only clear linkage that has been documented by our 
research is an effect that has originated outside of the rain forest and impinged upon the 
forest. We are yet to document effects of the rain forest upon the global atmosphere. 
We believe that we are in a better position now than before undertaking field 
measurements to speculate upon the potential effect of changes in the rain forest upon 
the global atmosphere. Further research is necessary before we can move from 
speculation to well-founded conclusions about the links and consequences between the 
rain forest and global system. 
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Abstract 


Meso- to synoptic-scale lines which form along the northern 
coast of South America as seabreeze-induced instability lines and 
propagate through the Amazon Basin are investigated using data 
collected during the April-May 1987 Amazon Boundary Layer 
Experiment (ABLE 2B) . 

These systems, termed "Amazon Coastal Squall Lines" (ACSL) , 
have been noted by others, but details of the structure and 
evolution of the ACSL are limited. The present paper uses 
analyses of GOES satellite, upper air rawinsonde and surface 
Portable Automated Mesonet (PAM) data to describe the structure, 
dynamics and life-cycle evolution of the ACSL. Twelve ACSL were 
sampled during ABLE 2B and three examples are discussed in 
detail . 

The ACSL are discontinuous lines of organized mesoscale 
cloud clusters which propagate across the Central Amazon Basin at 
speeds of 50-60 km h' 1 . The ACSL undergo six possible life cycle 
stages: coastal genesis, intensification, maturity, weakening, 

reintensification and dissipation. Analysis also indicates that 
the ACSL are composed of three distinct components: pre-storm 

cumulus (PSC) , leading edge convection (LEC) and multiple 
precipitating cloud layers in the trailing stratiform region 

(TSR) . 

Divergence and vertical velocity calculations indicate deep 
vertical ascent in the LEC and a region of mid-level convergence 
(w 500 mb) in the TSR. The latter is associated with a weak 
updraft above and an unsaturated downdraft below. The TSR 
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vertical motions are an order of magnitude smaller than in the 
LEC. 

Substantial shear in the low-level inflow occurred in all 3 
case studies and, as suggested by model simulations, may play an 
important role in the longevity (24-48 h) of the ACSL. Pre- and 
post-ACSL profiles of 0, demonstrate that the ACSL stabilize the 
troposphere in their wake and erode a pre-storm minimum ©-• It 
is hypothesized that the removal of this minimum is accomplished 
not only by direct mixing via vertical motions in the LEC ("hot 
towers") , but also through detrainment in the multiple-layered 
TSR. Part I describes the structure and kinematics of the Amazon 
Coastal Squall Lines (ACSL) while Part II deals with the heat and 
moisture transports of these systems. 
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1. introduction 

Among the most important rain producers in the rain forests 
of the Central Amazon Basin (CAB) of Brazil are meso- to 
synoptic-scale squall lines of coastal origin. These systems, 
which we will term the Amazon Coastal Squall Lines (ACSL) , form 
along the northeastern coast of South America as a result of sea 
breeze-induced convection. The ACSL often propagate away from 
the coast, reaching as far inland as Manaus in the CAB and, in 
some instances, reaching the Andes 48 h after genesis (Kousky, 
1980; Molion, 1987). 

Both Kousky (1980) and Molion (1987) have presented 
sequences of satellite images which illustrate the coastal 
genesis and inland propagation of a squall line in the Amazon 
basin. Molion (1987) found that these lines of instability, or 
squall lines, may be over 1000 km long and propagate across the 
Amazon basin with average speeds of 10° longitude per day or 45- 
55 km h' 1 . Kousky and Molion (1981) and Molion (1987) have 
suggested that the diurnal cycle and propagation characteristics 
associated with these squall lines may explain the relatively 
smaller precipitation values in the lower Amazon and the 
precipitation maximum in the central and western Amazon. 

More recently, Greco et al. (1990) used data collected 
during the wet season Amazon Boundary Layer Experiment (ABLE 2B) 
in April-May 1987 to investigate the ACSL, termed Coastal 
Occurring Systems (COS) in their work. Greco et al. (1990) found 
that these COS are generally 1000-2000 km long, but may attain a 
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length of 3500 km. In addition, they estimated that the COS 
propagated across the Basin at speeds of 50-60 km h l , and were 
observed to last 24-48 h. Similar results were reported by Cohen 
et al. (1989) in a study of "propagating squall lines" which 
utilized hourly satellite imagery data between 1979 and 1986. An 
important finding of Greco et al. (1990) was that 12 COS occurred 
during the ABLE 2B experiment (1 April-15 May 1987) and were 
responsible for 40% of the total measured rainfall. Cohen et al. 
(1989) found that these types of squall lines are most common in 

the CAB between April and August. 

The studies cited above provide only a qualitative 
description of the ACSL. Detailed investigations of the 
structure, flow field and energetics of the ACSL have been very 
limited (Scala et al., 1990) due mainly to the lack of data 
collected in the relatively harsh environment of the Amazon ram 
forest. The main goal of this paper is to use data collected 
from satellites and a mesoscale network of surface and upper air 
stations employed in the CAB during ABLE 2B to investigate the 
structure and kinematics of the ACSL on the meso- and convective 
scales. An additional goal is the comparison of the ACSL with 
other tropical and non-tropical squall lines. Although satellite 
imagery and preliminary analyses suggest similarities with the 
West African (Chalon et al., 1988; Chong et al., 1987), tropical 
Atlantic (Barnes and Sieckman, 1984), Venezuelan (Betts et al., 
1976) , and even the United States mid-latitude (Ogura and Liou, 
1980) squall lines, a more quantitative assessment is needed and 
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will be presented in this work. 

This paper, which is the first of a two part study on the 
ACSL , will document the life cycle and cloud components (section 

3) , and both the thermodynamic and kinematic structure (section 

4) of the ACSL as they propagate through the CAB. Important 
similarities and differences between the structure of the ACSL 
and other tropical and non-tropical squall lines will be 
discussed in Sections 4 and 5. Section 5 will also detail the 
important effects of the ACSL on the. larger scale environments of 
the CAB. Part II of this study will deal with the heat 
transports associated with the ACSL in the Amazon Basin. In 
particular, the focus will be upon the contributions of the ACSL 
heat transports to the total heat export from the tropics, and 
the relative importance of the heat transport produced by the 
stratiform regions of the ACSL. 

2. Data 

The National Aeronautics and Space Administration (NASA) , in 
cooperation with the Brazilian Instituto Nacional de Pesquisas 
Espaciases (INPE) and Instituto Nacional de Pesquisas da Amazonia 
(INPA) , conducted the wet season Amazon Boundary Layer Experiment 
(ABLE 2B) during 1 April-15 May 1987. A detailed description of 
the ABLE 2B can be found in Harriss et al. (1990) and Garstang et 
al. (1990). ABLE 2B, as part of NASA's Global Tropospheric 
Experiment (GTE) , was primarily an atmospheric chemistry 
experiment. However, ABLE 2B also included an extensive 
meteorological component which enabled the study of atmospheric 
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motions ranging from the turbulent to the planetary scale. 

The ABLE 2B meteorological networks provided the focus for 
meso- through synopt ic-scale sampling which was necessary for the 
study of convective systems over the vast, relatively 
inaccessible Amazon rain forests. The measurement systems 
operated during ABLE 2B included both ground and space-based 
platforms . 

Rawinsondes were released nearly simultaneously every 3 h at 
the three corner stations of a 1000 km 2 mesoscale triangle 
network located in the CAB (Fig. 1). Procedures for reducing and 
adjusting these data are described below. After data reduction, 
the soundings provided horizontal wind velocity, temperature, and 
mixing ratio from the surface to 100 mb. 

Portable Automated Mesonet (PAM) stations, the Geostationary 
Observational Environmental Satellite (GOES) platform and a j-cm 
radar acquisition system continuously monitored the atmosphere in 
the network triangle and the CAB. These measurements provided 
the data used to define and partition the ACSL systems into 
distinct cloud components. Four single level PAM stations (Fig. 
1), located 5 m above the rain forest canopy on 40 m towers, 
measured horizontal wind velocity, temperature, humidity, 
pressure and precipitation within the mesoscale network. The 
PAM measurements, obtained at a frequency of 1 Hz and stored as 1 
min averages, were transmitted via satellite to the National 
Center for Atmospheric Research (NCAR) for processing and 


storage . 
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The GOES data consisted of visible and infrared (IR) images 
obtained via a Digital Weather Image Processing System (DWIPS) 
operated at the ABLE 2B operations center at Eduardo Gomes 
Airport (EGA) in Manaus, Brazil. The GOES images were used to 
analyze cloud conditions over the mesoscale triangle and to 
document the movement of synoptic disturbances. An ACSL life- 
cycle taxonomy is derived from images obtained over the entire 
experiment . 

A 3 cm radar located at EGA provided Constant Altitude 
Planned Position Indicator (CAPPI) scans of the mesoscale 
triangle at 5 min intervals. Signal attenuation severely 
degraded the CAPPI data during heavy precipitation, thereby 
limiting the radar's usefulness for making accurate estimates of 

rainfall rates. 

The mesoscale rawinsonde data used in this study were 
checked against climatology and rawinsonde profiles measured 
coincidently by a 6-station basin-scale network (Fig. 1) which 
was also part of the ABLE 2B experiment. This comparison 
revealed pronounced pressure offset problems in the mesoscale 
triangle sounding which were ultimately traced to bad instrument 
calibrations by the rawinsonde manufacturer . 

A least squares procedure based on the mean basin scale 
temperature profile was used to develop a procedure for adjusting 
and editing the mesoscale soundings. After the pressure offset 
was determined, temperatures were reassigned to adjusted pressure 
levels. Geopotential, mixing ratio, and wind velocity were then 
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recalculated at 25 mb intervals from the surface to 100 mb. The 
wind calculations employed a smoothed azimuth, elevation, and 
balloon ascent rate derived from a three-point moving average 
(Davis, 1986) to reduce noise introduced by antenna errors. 

Some soundings launched during thunderstorms reflected data 
losses due to icing problems. Usually these soundings terminated 
near the freezing level. These weather-related problems 
prevented budget calculations for several strong ACSL which 
passed through the mesoscale network. The cases selected for 
study either had no problems or contained data deficiencies which 
were easily handled by the above procedures. 

Divergence for the mesoscale triangle was computed from the 
line integral of the wind components normal to each leg (Gamache 
and Houze, 1982) . Divergence was calculated from 1000 to 100 mb 
in 25 mb intervals. The vertical velocity at each pressure level 
was determined first from the vertically integrated divergence, 
and then adjusted by O'Brien's (1970) method which forces the 
vertical velocity to zero at the surface and at the highest level 

(100 mb) . 

This procedure gives results that compare well with 
physically realistic estimates of vertical motion from adiabatic 
and isentropic methods (Haltiner and Martin, 1957) . 0 Brien s 

method assumes vertical velocity errors result from linear 
divergence errors which increase uniformly from a surface minimum 
to a maximum at the highest level reached by the sonde. Gamache 
and Houze (1982) obtained acceptable results by applying these 
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procedures to tropical squall lines. 

3. Results 

a . Definition of the ACSL L ife Cycle 

Analysis of both GOES imagery and PAM surface data have 
shown that 12 ACSL were detected in the CAB during the ABLE 2B 
wet season (Greco et al., 1990). These large and generally 
linear meso- to synoptic scale systems form along the northern 
coast of South America in mid-afternoon and move across the 
Amazon Basin, usually from NE to SW, at an- average speed of 50-60 
km h' 1 . The ACSL often last between 24 and 48 h, reaching the 
ABLE 2B experiment region in the central Amazon Basin 20—24 h 
after formation along the coast. Figure 2, taken from Greco et 
al. (1990) , shows the composite hourly cloud cover percentages 
from the GOES-W brightness temperatures less than -33 °C for the 
12 ACSL days. The development of the squall or cloud line along 
the coast by 2200 UTC (1800 LST ) is clearly shown. The 
existence, at the same time, of another cloud line in the CAB 
should be noted. This is the previous day's ACSL. The 
propagation of the cloud line or ACSL from near the coast at 2000 
UTC to the CAB 18 h later (1600 UTC) is graphically shown. 

At maturation, ACSL are usually 1000-2000 km long and 100- 
300 km wide, but can reach a total length of over 3000 km and are 
often a dominant synoptic scale cloud feature over the third of 
the tropics appearing on GOES full-disk images (Molion, 1987). 
This is illustrated by Fig. 3 which shows the ACSL of 26 April 
1987 as it approaches the mesoscale network. This ACSL was 
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approximately 3500 km long at the time of strongest intensity, 
although only 1000 km may have contained active convection. The 
ACSL is not a continuous line of single cells, but rather a 
highly discontinuous arc or line of discrete clusters of cells. 
Although organized on the mesoscale, the clusters vary in their 
dimensions within the line. Each mesoscale cluster is composed 
of deep cumulonimbus along the leading edge with a shield of 
anvil cloud that surrounds the convective cells and extends far 
to the rear. 

Greco et al. (1990) showed that these massive, long-lived 
convective systems accounted for at least 40% of the total 
rainfall in the mesoscale network during ABLE 2B. Their 
importance to the water budget of the Amazon Basin makes it 
necessary to undertake a detailed study of their life cycle and 
the evolution of the flow field and kinematic structure. 

The GOES visible and IR images collected during ABLE 2B 
indicate that patterns, sizes, intensities of clouds and rainfall 
can change dramatically during the lifetime of a ACSL. The GOES 
images from April and May 1937 have been studied to define a ACSL 
life cycle of six possible stages. These stages are: 

1) COASTAL GENESIS 

During coastal genesis, the first thunderstorm cells form 
during the afternoon in the seabreeze convergence zone along the 
northern coast of South America. The number and size of the 
cells increases -with time, followed by selective merging into 
thunderstorm clusters with shared anvils. 
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2) INTENSIFICATION 

During late afternoon and evening, merging continues as 
small clusters of thunderstorms merge to form larger clusters, 
generally oriented from northwest to southeast in a quasi-linear 
pattern parallel to the coastline of NE Brazil. Figure 4a 
illustrates the intensification of clustering of the sea breeze 
forced convection on the evening of 25 April 1987. The strongest 
clusters containing intensifying convective elements and anvils 
are resolved by segment number 7 (black) on GOES images using the 
Mb enhancement technique (Carlson, 1981) . The strong clusters 
have sharp boundaries and cloud-top temperatures colder than 
-59 . 2 °C , indicating a height of approximately 13-14 km. 

3) MATURATION 

During this stage, the entire system achieves maximum length 
and width and propagates at average speeds of 45-55 km h‘ l . 
Individual clusters also reach their largest size and 
strongest intensity, as seen on the enhanced IR images. A mature 
ACSL may cover a total length of 2000-3000 km with 30-40% of the 
leading edge consisting of active deep convection. The width of 
a mature ACSL is typically 200-300 km. Patches of white (Mb 
segment 9) in the mesoscale clusters indicate cumulonimbi with 
cloud top temperatures less than -80.2°C and heights of 16-17 km. 
The merging clouds between individual clusters give the ACSL 
distinctive synoptic scale features on the GOES full disk images. 
This is illustrated in Fig. 4b, which shows an ACSL that formed 
on 25 April and had already reached maturity over the CAB by 
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early afternoon of 26 April. The formation of a possible 
new ACSL is also evident along the coast at this time. 

4) WEAKENING 

The clusters diminish in size and their outer definition 
becomes ragged on the Mb enhanced images. Areas previously shown 
in black on the Mb images change to dark and light gray (segment 
numbers 6 and 5, respectively) , indicating warmer cloud-top 
temperatures and weaker thunderstorms. The width of the ACSL 
decreases to 100-200 km. 

5) REGENERATION 

During the following afternoon, ACSL regeneration resulting 
from increased daytime surface heating, may follow weakening 
which occurred at night. During regeneration, an enlargement of 
the black areas (segment 7) with a sharp outer definition on the 
Mb images indicates thunderstorm clusters that are growing and 
intensifying. The system's synoptic-scale features become better 
defined on the full disk IR images. Figure 4c shows the ACSL of 
25-26 April as it regenerates in the western Amazon Basin at 0300 

UTC 27 April. 

6) DISSIPATION 

After a period of weakening, the individual clusters become 
more ragged and ultimately disappear. The system loses its 
synoptic scale definition as shown on the full disk images. 

The above definition of the ACSL life cycle is similar to 
taxonomies developed for tropical squall line-type mesoscale 
convective systems (e.g., Houze, 1977; Frank, 1978; Leary and 
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Houze, 1979) . This present taxonomy, however, defines the life 
cycle of a meso- to synoptic-scale convective disturbance. 
Utilizing this life cycle taxonomy, we have classified three ACSL 
as they passed through the mesoscale network; the ACSL on 26 
April 1987 was classified as weakening, the one on 1 May 1987 as 
mature, and the one on 6 May 1987 as a regenerating system, 
b . ACSL Cloud and Rainfall C omponents 

The ACSL have a convective and stratiform structure similar 
to the design of Gamache and Houze (1982) and satisfy Zipser's 
(1977) criteria for a tropical squall system. Figure 5, taken 
from Garstang et al. (1990), shows a schematic representation of 
a typical mature ACSL. The three primary cloud components 
comprising ACSL are: pre-storm cumulus (PSC) ahead of the squall 

front, leading edge convection (LEC) , and widespread layers of 
anvil clouds in the trailing stratiform region (TSR) . 

The three cloud components of each ACSL were discerned 
through subjective analysis of both half-hourly GOES images with 
a 2 km resolution and 5 min CAPPI radar scans also with a 2 km 
resolution. Most importantly, however, was an analysis of the 
high frequency PAM data which was used to delineate the three 
ACSL components. The 1 min averaged values of T, q, © ; and wind 
speed clearly identified the passing of the squall front. 
Rainfall rates measured by the tipping bucket raingauge mounted 
on the PAM platform were employed to determine the transition 
from deep convection stage to anvil stage using .25 mm min 
cutoff from convective to stratiform rain. The .25 mm min 
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cutoff was necessitated by the .25 mm resolution of the tipping 
bucket raingauge. 

Figure 6 shows the time series of PAM observations recorded 
at Embrapa on 26 April 1987. The arrival of the gust front is 
clearly noted between 1630 and 1700 UTC. The marked response of 
selected surface variables (as measured by the PAM network) to 
the passage of the gust front at Ducke, Embrapa, Carapana and ZF 
1 on 26 April 1987 and 6 May 1987 are summarized in Table 1. 
Results for 1 May 1987 are not . shown because lightning disrupted 
the PAM network during ACSL passage. There are large differences 
among the four PAM stations in rainfall rate and changes of T, q, 

P and 9, with the gust front passage, especially on 6 May 1987. 
These differences reflect the spatially discontinuous nature of 
the ACSL. 

4. Squall Line Kinematics and Dynamics 

In this section, analysis of the mesoscale network PAM and 
rawinsonde data is used to reveal the cross-sectional (normal to 
direction of storm propagation) structure and characteristics of 
the 1 May 1987 ACSL. We have chosen to emphasize details of the 
1 Mav case because it was the only one of the three for which 
rawinsonde coverage existed for both the LEC and TSR components . 
Observations taken from the mesoscale triangle, however, did 
capture the TSR components of the 26 April and 6 May systems, and 
we compare the general kinematic structure of this region for the 
three squall systems. Model simulations of the 6 May and 26 
April squall lines have been presented by Scala et al. (1990, 
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1991 ) . 

The 1 May ACSL formed during the mid-afternoon of 30 April 
1987 as a line of numerous small convective clusters located 
along the northeast coast of Brazil. By late afternoon, the 
clusters had grown rapidly and merged to form a synoptic arc- 
shaped ACSL. By early evening, the arc acquired a more linear 
appearance as it propagated toward the southwest at 33 km h 
(Fig. 7 ) ■ A group of smaller mesoscale clusters formed 250 km 
ahead of the ACSL leading edge. These . clusters produced heavy 
rain in the northern part of the mesoscale triangle six hours 
prior to the arrival of the ACSL at 1200 UTC, 1 May. 

Mass divergence was calculated over the mesoscale network 
volume using the near-simultaneous rawinsonde launches at the 
triangle vertices (Ducke, Embrapa, Carapaha) and by employing a 
simple line integral method. Vertical profiles of the divergence 
values calculated for the LEC (1200 UTC) and TSR (1800 UTC) 
components of the 1 May ACSL are shown in Fig. 8. The LEC 
profile is characterized by a pronounced vertical couplet of 
convergence (-1.6 x 10' 4 s' 1 ) between 600-7 00 mb and divergence 
(1.8 x lO" s' ) located between 300-500 mb within the deep 
convective region. Vertical profiles of vertical p-velocity for 
the LEC and TSR components are presented in Fig. 9. The 
divergence couplet is associated with deep upward motion in the 
LEC with a peak of 40 /libs' 1 near 550 mb. Figure 10 shows the pre- 
storm (1200 UTC) and post-storm (1300 UTC) composite vertical 
profiles of temperature and dewpoint for the triangle. Note that 
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the deep convective ascent coincides with saturated conditions to 
a level of 500 mb. It is likely that the 1200 UTC sounding 
ascended through and directly sampled a portion of the deep 
convective towers of the squall disturbance. 

Six hours following the storm passage, the divergence 
profile of the TSR is relatively weak throughout the depth of the 
troposphere and includes a divergent layer that extends from 800 
mb to the surface. Peak downward motion is indicated at 800 mb 
in the TSR, accompanied by- progressive warming and drying of the 
layer (see Figs. 9 and 10). Higher in the atmosphere, weak 
ascent (5-10 Mbs') is found between 600 mb and the tropopause. 

The peak in the upward motion coincides with the upper portion or 
a mid-level saturated (cloud) layer seen at 1800 UTC in the 
thermodynamic profiles. While the 1800 UTC sounding indicated 
that the boundary layer was recovering from the convective 
disturbance (note the surface layer warming and moistening 
between 1200 and 1800 UTC) , the formation of a warm, dry 
subsidence-induced inversion in the low levels combined with a 
distinct region of warming above 500 mb was sufficient to 
suppress deep convective activity in the wake of the squall line. 

The vertical p-velocity values are consistent with the 
findings of Gamache and Houze (1982) who have shown that upward 
vertical motion is an order of magnitude larger in the convective 
region than in the stratiform region. The updraft structure 
comprised of a concentrated core of rising motion in the leading 
convective elements and weaker ascent in the trailing stratiform 
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cloud region is a feature that characterizes many mesoscale 
convective systems in the tropics, including systems observed 
during COPT81 (Chalon et al., 1988), AMEX (Frank and McBride, 
1989), GATE (Gamache and Houze, 1982), WMONEX (Johnson and 
Priegnitz, 1981); and in the midlatitudes (Gallus and Johnson, 

1991) . 

Figure 11 presents a comparison of vertical p-velocity 
calculations for each of the three TSR components on 26 April, 1 
May and 6 May. In general., weak ascent occupies the upper 
portion of the troposphere (above 500 mb) in each of the three 
TSR components, with descending motion in the lower layers. The 
weakest motions are those associated with the 6 May regenerating 
case. It should be noted that the upward motion in the weakening 
26 April system occupies nearly the entire tropopause (except for 
a shallow subsidence layer below 850 mb) , and the magnitude of 
the motions in this deep layer are nearly the same as that for 
the mature system. These findings may suggest that ascent in the 
trailing anvil plays a more dominant role in the late stages of 
the ACSL life cycle. However, since our analysis compares the 
life cycle stages of individual systems on different days, the 
observed differences in velocity structure might also reflect 
variability in the large scale tropospheric conditions which 
influence the intensity and organization of the sguall line 
circulations . 

The pre-storm line-normal component of the relative flow for 
each of the three ACSL was also calculated. These calculations 
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utilize a storm coordinate system defined by an x' axis oriented 
orthogonal to the orientation of the leading edge of each squall 
line (Fig. 12) . The angle alpha between the storm coordinate x' 
axis and the x axis in the standard meteorological reference 
frame is used to derive a u— wind component normal to the squall 
orientation (u'-wind component) . The u' wind component was taken 
from pre-storm triangle composite profiles of wind speed. These 
profiles are shown in Fig. 13. The storm forward speed was then 
subtracted from this u'-wind component to obtain the line-normal 
relative flow. For the three cases of 26 April, 1 May and 6 May 
we have defined alpha as 50°, 60° and 45°, and used storm speeds 

of 6 m s' 1 , 9 m s’ 1 and 10 m s' 1 , respectively. 

The pre-storm relative flow profiles for all three cases are 
presented in Fig. 14. All three cases show a pronounced shear m 
the lowest levels. This rear to front low level inflow is mainly 
a result of the easterly jet that characterizes the lower wet 
season troposphere over the Central Amazon Basin (located at 2.5 
km in Fig. 13). A similar low-level jet structure, the African 
easterly jet, is associated with West African squall lines 
investigated during COPT 81 (Chong et al., 1987). The observed 
relative flow profiles from COPT 81 are very similar to those 

identified for Amazon systems. 

The low-level shear associated with the three ACSL also 
suggests a mechanism for their observed longevity. Rotunno et 
al. (1988) have demonstrated that the vigorous, sustained lifting 
along the squall leading edge is maintained by a balance between 
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counter-rotating circulations generated by the shear in the 
outflow cold pools and the opposing low-level relative inflow. 
Thorpe et al . ( 1982 ) have proposed that strong low— level shear 

relative to the storm prevents the outflow boundary from 
propagating away from the storm leading edge and thus cutting off 
the storm's available supply of moist boundary layer air. The 2- 
D version of the Goddard Cumulus Ensemble (GCE) model (Tao et 
al . , 1991) has been used to simulate the Amazon systems (Scala et 

al. , 1990, 1991) . 

The results of efforts to simulate the 26 April and 6 May 
convective events also indicated the potential for convective 
demise when the updraft and downdraft interface sloped too far 
upshear. This action initiated a propagating outflow boundary 
that terminated the influx of high Q e air to the developing 
cells. A convective balance between the low-level shear and 
storm— induced outflows is likely reguired for the long-lived 
storms observed over the Amazon basin. 

The distinct minimum in relative flow seen at 2.5 km for the 
three Amazon systems is a direct consequence of the easterly jet 
and, in at least two cases (26 April and 6 May), indicates the 
presence of weak inflow into the rear of the convective region. 
This inflow can provide a source of momentum and potentially low 
©. air for convective— sea le downdrafts. Low 0 e air also enters 
the front of the three squall cases at mid-levels between 4 and 8 
km. The trajectory of this air within the storm system is 
uncertain, although model results (Scala et al., 1990, 1991, 
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1992 ) suggest parcels can be entrained in the updraft cores or 
pass through them to the rear where they may be incorporated in 
convective scale downdrafts. Houze ( 1977 ) suggests that this 
type of inflow may in fact pass around updraft cores and 
contribute to the production of convective downdrafts on the 
trailing sides of the system. The mid-level inflow may also 
transport hydrometeors that are generated in the convective 
towers into the trailing portions of the system. These 
horizontal transports have' been found to account for a 
substantial portion of the condensate that comprises the trailing 
stratiform clouds of squall lines (Tao et al., 1991; Chong and 
Hauser, 1989) . 

Differences in the strength of the mid-level inflow between 
the three Amazon systems may be a consequence of the variation in 
the larger-scale tropospheric flow regime, either as variations 
of the day-to-day flow, or from a modification of the basic state 
flow produced by passage of a large squall line during the 
previous day. 

The relative flow profiles also demonstrate somewhat 
different characteristics in the uppermost levels. Studies with 
the two-dimensional cloud model (Scala et al., 1991, 1992) show 

that modification of updraft cores by the environmental wind 
field produces one of two types of anvil circulation, either an 
overturning or a trailing type of structure. On the one extreme, 
the environmental flow found at 15-16 km on 1 May promotes 
development of an overturning anvil, while the strong upper level 
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westerlies on 6 May favor development of a trailing anvil 
structure . 

As the different features observed in the relative flow 
fields for the 26 April, 1 May and 6 May squall systems 
demonstrate, the gross structure and behavior of convective 
systems in the ACSL classification may be similar, but variations 
in the finer aspects of the storm flow do occur. These 
variations reflect both differences in the large scale flow 
regime and sampling of the particular stage of the storm life 

cycle . 

5. Effects on the Larger Scale Environment 

Pre- and post-storm vertical profiles of 0 e for the 26 
April, 1 May and 6 May ACSL are shown in Fig. 15. The pre-storm 
profiles on 26 April, 1 May and 6 May were taken at 1-2 hours 
before the passage of the ACSL. The three pre-storm profiles all 
display a minimum in © c of about 339°K. Differences in the 
overall shape of the pre-storm profiles in the lower troposphere 
are due primarily to the time of day and the amount of associated 
surface heating, although the isothermal structure of the lower 
portion of the 1 May profile may reflect the actual ascent of the 
rawinsonde through the core of deep convective towers. 

Comparison of the pre-storm and post-storm profiles reveals 
that the Amazon systems have a pronounced effect on their 
tropospheric environment. These systems effectively reduce the 
mid-level pre-storm minimum in 0.. This is accomplished by the 
upward transport of surface air containing high 0 S values and the 
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release of latent heat in mid-levels, and the compensating 
downward transport of cool, dry mid-tropospheric air. Evidence 
for mid-level warming is seen in the soundings shown in Fig. 10. 
The troposphere in the wake of the 1 May system has been warmed 
nearly 2°K between 525 and 250 mb. An accompanying moistening of 
nearly 1 g kg' 1 has also occurred. Compensating downward 
transport of dry (relatively low 0 e ) air is also noted through 
the formation of a pronounced dry inversion between 700 and 950 
mb. This stabilization of the vertical column in the wake of the 
system is noted in the profiles for each of the three Amazon 

cases . 

The profiles of 9. suggest that the troposphere is well 
mixed by the ACSL.' This finding, also concluded by Scala et al. 
(1990) , implies that the Riehl and Simpson (1979) "hot tower 
requirement for vertical transport may not be operating within 
these systems over the Amazon Basin. Transport through protected 
hot tower conduits, if operating, would preserve the 9 e minimum 
in the wake of the squall line. Instead, erosion of the 0.. 
minimum and low level decreases in 0, imply that vertical motions 
in the ACSL directly overturn the troposphere. Vertical 
transports of heat under these circumstances can be accomplished 
by mean upward and downward circulations ranging in scale from 
that of the ACSL to an integrated Hadley transport. The cores of 
the deep convective elements may in fact "leak" or actively 
detrain latent and sensible heats into the surrounding 

Evidence for this process was seen during aircraft 


environment . 
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flights which revealed multiple layers of cloud that formed in 
the rear of the deep convective towers. The entrainment of 
relatively dry mid-tropospheric air into the updraft cores and 
subsequent dilution and detrainment to the rear of the system 
provides a possible mechanism for the "leaking" process. 

The 0„ profiles for 1 May also reveal that upper- 
tropospheric cooling occurred in the wake of this system. This 
is reflected in the post-storm sounding by a temperature 
reduction of over 3°K above 175 mb. As indicated by the profiles 
of ©,, the combination of mid-tropospheric warming and this 
observed upper-level cooling renders the middle and upper 
troposphere less stable in the wake of the storm. The origin of 
the cooling is uncertain (perhaps due to ice sublimation) , but we 
note that cooler conditions occurred aloft in the tropospheric 
environment ahead of the squall line. The 1 May 0, profile is up 
to 2°K cooler than the 6 May profile and 4°K cooler than the z6 
April profile. The reduction in stability shown in the 0 4 
profiles may help sustain the observed upward motions in the TSR. 
5. Summary and Conclusions 

In this paper , we have examined the structure and dynamics 
of a distinct type of synoptic-scale convective system over the 
equatorial continental tropics, the Amazon Coastal Squall Line 
(ACSL) . ACSL are discontinuous lines of organized mesoscale 
cloud clusters that form along the northern coast of South 
America as seabreeze- induced instability lines. The lines then 
propagate across the Central Amazon Basin at speeds of 50 60 km 


h' 1 with some maintaining enough coherence to arrive at the 
westernmost boundaries of the Amazon Basin 24-48 h later. 

Analysis of satellite and rawinsonde data have shown that 
the ACSL undergo six possible stages in their life cycle, 
genesis, intensification, maturity, weakening, reintensif ication 
and dissipation. Detailed analyses of surface PAM, satellite and 
aircraft observations also reveal that ACSL are composed of three 
distinct components: pre-storm cumulus (PSC) , leading edge 

convection (LEC) , and multiple, precipitating layers of cloud in 
the trailing stratiform region (TSR) . 

Calculations of divergence and vertical velocity for a 
mature ACSL indicate deep vertical ascent in the LEC and a region 
of mid-level convergence in the TSR. This convergence is 
associated with a weak updraft above 500 mb and unsaturated 
downdraft below. The vertical motions in the TSR are an order of 
magnitude smaller than those in the LEC. These findings are 
similar to the structure and dynamics of squall lines and other 
types of mesoscale convective system in the tropics and 
midlatitudes (Houze, 1989) . 

Substantial shear in the low-level inflow was found to occur 
with the 26 April, 1 May and 6 May squall lines. The shear is a 
consequence of an easterly wind maximum located at 2 . 5 km, and 
simulations using a two-dimensional cloud model (Scala et al., 
1990, 1992) suggest that the shear may play an important role in 

the longevity of the storm system. Mid-level inflow of air into 
the leading convective elements was found to vary widely between 
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the three systems, and may act to dilute the updraft cores as 
well as advect appreciable quantities of condensate into the 
trailing portions of the squall system. 

pj-of j_ ]_es of pre- storm and post - storm Q c demonstrate the 
effect of ACSL on their larger-scale tropospheric environment. 

The ACSL generally stabilize the troposphere in their wake and 
erode the pre-storm minimum in Q c . We hypothesize that the 
removal of this minimum is accomplished through direct mixing via 
vertical motions in the storm, and also through, detrainment in 
multiple layers in the wake of the squall system. The evidence 
supports the contention that energy transport through protected 
"hot towers" is not the preferred mechanism that accomplishes the 
required vertical transports over the Amazon Basin during the wet 
season. Energy budgets and transports associated with the ACSL 
will be the main focus for Part II of our study. 
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Fig. 1. 


Fig. 2. 


Fig. 3. 


Fig. 4. 


Fig. 5. 


Fig. 6. 


List of Figures 

Location and design of the meteorological network used 
from 1 April-15 May 1987 during ABLE 2B. 

Composite hourly cloud cover percentage for GOES-W 
brightness temperature < 33°C for 12 ACSL days during 
ABLE 2B. Times shown are (a) 2200 UTC, (b) 0400 UTC, 

(c) 1000 UTC and (d) 1600 UTC. 

GOES full-disk infrared image of the 26 April 1987 
ACSL as it crossed the Central Amazon Basin. White 
circle in center of image depicts the location of the 
instrumented mesoscale triangle. 

Sequence of GOES full-disk images showing life cycle 
stages of the 25-27 April 1987 ACSL. Times shown are 
(a) 2100 UTC, 25 April, (b) 1800 UTC, 26 April and (c) 

0300 UTC, 27 April. The + symbol denotes the location 
of Manaus, Brazil. 

Schematic representation of cloud and weather features 
associated with a mature ACSL during ABLE 2B. The 
features are based on observations made from the 
surface, airborne platform and satellite imagery. 
Airflow pattern is inferred from relative flow 
and vertical velocity calculations in the triangle 
volume . 

Time series of PAM measured surface variables during 
passage of the 26 April 1987 ACSL. One-minute averages 
are shown for (a) rainfall rate in mm min* 1 , (b) peak 
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wind gust in m s' 1 , (c) u-wind component in m s' 1 , (d) v 

wind component in m s 1 , (e) piressuire in hPa, (f) 

temperature in degrees (C) , (g) mixing ratio in g kg 


Fig. 7. 


Fig. 3. 


Fig . 9 . 


Fig. 10. 


Fig. 11. 


and 0. in K. 

GOES Mb-enhanced infrared image of the 1 May 1987 
ACSL at 0230 UTC as it propagated westward across the 
mesoscale triangle (marked by the +) . The letters G 
and B denote the location of, respectively, Georgetown, 
Guyana and Belem, Brazil. 

Vertical profiles of mass divergence on 1 May 1987. 
Solid curve depicts the sguall line LEC component a<_ 
1200 UTC. Dashed curve is the TSR component at 1300 
UTC. Values are in units of 10° s' 1 . 

Vertical profiles of vertical p— velocity for the 1 
May 1987 ACSL. Solid curve depicts the squall line 
LEC component at 1200 UTC; dashed curve is the TSR 
component at 1800 UTC. Values are in units of M b s' 1 . 
Vertical profiles of sounding-derived temperature and 
dewpoint for the 1 May 1987 ACSL in standard skew-T/ 
log P format. Mesoscale triangle composite soundings 
are used for pre-squall (1200 UTC, solid curve) and 
post-squall (1800 UTC, dashed curve) observations. 
Vertical profiles of vertical p-velocity for the TSR 
of three squall lines taken at 1300 UTC. TSR diver- 
gence values correspond to 26 April (solid curve) , 

1 May (dashed curve) and 6 May (dash-dot curve) 19o7 
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Fig. 12. 


Fig. 13. 


Fig. 14. 


Fig. 15. 


squall lines. 

Coordinate system used to derive the relative flow 
line— normal component for the 1 May 1987 ACSL. 
Profiles of rawinsonde-derived wind velocity for 
three ACSL pre-storm environments: 26 April 1987 

(1500 UTC) , 1 May 1987 (1200 UTC) and 6 May 1987 

(1500 UTC). Wind direction is according to standard 
compass . A half-barb equals 2.5ms; a full barb 
equals 5 m s' 1 . 

Horizontal storm relative flow calculated for the 
pre-storm components of the 26 April, 1 May and 6 May 
ACSL. Wind speed is shown in m s' 1 . 

Profiles of pre-storm and post-storm 0„ derived for 
the triangle composite soundings for three ACSL: 26 

April 1937 (1500 and 1800 UTC) , 1 May 1987 (1200 and 
1800 UTC), and 6 May 1987 (1500 and 1800 UTC). 


TABLE 1. 


values of temperature, humidity, © 


Summary of changes in surface 
and pressure observed at the 4 PAM locations with the passage of the squal 
line on (a) 26 April and (b) 6 May. Also presented are peak wind gusts, 

maximum rain rate and both total amount and duration of rainfall. 


a) 26 April 1987 

Ducke Carapana Embrapa 


AT 

-6 0 C 

-6. 5°C 

-7 . 0°C 

i 

cr» 

o 

o 

O 

Aq 

-3.7 g kg- 

1 -2.1 g kg 

1 -2.5 g kg' 1 

-2.6 g kg'* 

A©, 

-17.0 K 

-11.0 K 

-16.0 K 

— 

AP 

+0.7 hPa 

4-0.8 hPa 

+0.9 hPa 

+0.9 hPa 

Peak Gust 

13 m s* 1 

9 m s' 1 

12 m s' 1 

13 m s' 1 

Maximum Rain 





Rate 

0.25 mm min' 1 

1.5 mm min' 1 

1.7 mm min'* 

1 mm min' 1 

Total Rainfall 

1 mm 

3 4 mm 
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Abstract 


The role of the Amazon basin as a global center of action 
for deep convection is quantified in the context of the planet's 
heat balance. Mesoscale heat and moisture budgets are presented 
for three synoptic scale Amazon coastal squall lines (ACSL) over 
the central Amazon basin during the 1987 wet season. The budgets 
are obtained from volumetric analyses of mesoscale rawinsonde 
data from the second Amazon Boundary Layer Experiment (ABLE 2B) . 
Classification of cloud and rainfall components are based on data 
from Portable Automated Mesonet (PAM) towers, GOES images, and 
radar. Satellite estimates of the total active cloud area, 
coupled with the mesoscale budget calculations, are used to 
determine the upper level vertical heat transport for an entire 

ACSL as a function of life cycle. 

The instantaneous upper level heat transport in a mature 
ACSL equals nearly 20% of the theoretical export requirement for 
the equatorial trough zone. A little more than half of the 
system-wide transport occurs in the dynamically active anvil 
cloud. Amazon coastal squall lines are shown to exist over the 
basin simultaneously in various stages of formation, maturity and 
decay. When three stages are present, the combined instantaneous 
production of heat exceeds one third of the requirement to 
balance the global heat budget. 
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1. introduction 

Amazon Coastal Squall Lines (ACSL) , classified as "Coastal 
Occurring Systems (COS)" by Greco et al. (1990), "Propagating 
Squall Lines (PSL) " by Cohen et al. (1989), and identified as 
lines of instability by Kousky and Molion (1981), rank among the 
largest of tropical synoptic disturbances. These systems can 
reach a total length of 2000-3000 km, with over 1000 km being 
active deep convection. More importantly, Greco et al. (1990) 
found that these "COS" were responsible for at least 40% of the 
total precipitation which fell over the Central Amazon Basin 
(CAB) during the 1987 wet season months of April and May. 

In light of this important contribution to rainfall 
production, Part I of this study investigated the life cycle and 
kinematic structure of the ACSL over the CAB. As detailed in 
Part I, the ACSL develop during late morning-early afternoon 
along the northern coast of South America within an area of 
seabreeze-induced convection. They then intensify and propagate 
inland, possibly as mesoscale waves associated with trapeze 
instablity (Sun and Orlanski, 1981a, b) . At night, due to the 
cutoff of surface heating, the ACSL and its associated convection 
weaken. Many of the ACSL dissipate before they enter the CAB. 
However, a large number (12 during the 1 April-15 May 1987 Amazon 
Boundary Layer Experiment) of ACSL continue propagating through 
the night and regenerate the following day in the CAB due to 
increased surface heating. The average propagation speed of the 
ACSL before dissipation is 50-65 km h' 1 . At maturation, ACSL are 
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usually 1000-2000 km long and 100-300 km wide. These results 
agree well with the findings of Kousky (1980), Kousky and Molion 
(1981), Cohen et al. (1989) and Greco et al. (1990). 

Divergence, vertical velocity and storm-relative flow 
calculations presented in Part I suggest a primarily 2-D airflow 
structure which features deep upward motion in the leading edge 
convection (LEC) and a region of mid-level speed convergence in 
the trailing stratiform region (TSR) . Mesoscale ascent occurs in 
the TSR above 500 mb and subsidence- below. The vertical motions 
in the TSR are an order of magnitude smaller than those in the 
LEC. These findings agree well with the structure and dynamics 
of other squall lines found both in the tropics and midlatitudes 

(Houze , 1989 ) . 

Part II of this study evaluates the mesoscale heat and 
moisture transports associated with three ACSL sampled in the CAB 
during the 1987 Amazon Boundary Layer Experiment (ABLE 2B) . 

Special emphasis will be placed upon investigating the evolution 
of the heat and moisture transports as a function of both ACSL 
life cycle and cloud component. A secondary focus will be an 
evaluation of the relative importance of ACSL heat transports in 
the context of the total required global export of heat from the 
tropics. 

The Amazon Basin plays a recognized role as a global "center 
of action" for thunderstorms (Newell et al., 1972), with the 
storms often concentrated in synoptic disturbances (Kousky, 1979, 
1980; Molion, 1987; Greco et al., 1990). Until now, insufficient 
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rawinsonde measurements from the vast, relatively inaccessible 
Brazilian rain forests have precluded a definitive treatment of 
ACSL structure and energetics. This study uses a mesoscale 
network of surface and rawinsonde observations from the Apnl-May 
1987 ABLE 2B field program (Garstang et al. , 1990; Harriss et 

al., 1990) to assess ACSL energetics. 

While both Riehl and Malkus (1958) and Riehl and Simpson 
(1979) invoke the "hot" tower hypothesis to account for the 
needed upward transfer of heat from the lower to upper atmosphere 
in the tropics, they do not emphasize the concentration of the 
towers within any specific region of the tropics. They state 
that the required convective vertical transport of heat occurs in 
1500-5000 giant clouds active at one time and distributed in 
approximately 30 large disturbances throughout the tropics. 

The results presented in this paper suggest that organized 
meso- to synoptic-scale systems, such as the ACSL, which occur 
over a large equatorial, continental land mass contribute a 
significant fraction of the heat required to balance the heat 
budget of the planet. 

The findings reported in the following sections differ from 
the results of Riehl and Malkus (1958) and Riehl and Simpson 
(1979) in that we suggest that injections of heat in the 

equatorial trough can be: 

a) concentrated in centers of action 
land mass such as South America. 


over a continental 
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b) concentrated in organized meso- to synoptic-scale 
systems in which only half of the vertical heat 
transport is by the "hot tower" mechanism and half 
by mesoscale stratiform circulation. 

The ACSL are further shown to generate multiple stratiform cloud 
layers in the moist environment of the CAB. Mesoscale 
circulations associated with the stratiform cloud layers are 
shown to accomplish up to 50% of heat transported to the upper 
troposphere. Hot tower transport occurring in the deep 
convective towers accounts for the remaining 50%. 

Recent budget studies of tropical squall lines (Johnson, 

1984; Houze, 1989; Ackerman et al., 1988; Gallus and Johnson, 

1991) have also substantiated and clarified the important role 
played by anvil cloud layers. By partitioning their squall line 
analysis into convective and stratiform components, they show 
marked mesoscale ascent and vertical transports in the stratiform 
anvil. The energy needed for this mesoscale ascent has been 
attributed to the horizontal transport and injection of 
hydrometeors from the convective regions into the trailing 
stratiform layers (Rutledge and Houze, 1987). The total amount 
of heating will be shown to vary as a function of anvil area 
during the ACSL life cycle. 

2 . Data 

The wet season ABLE 2B was in operation over the CAB from 1 
April to 15 May 1987. A detailed description of the type, 
quality and frequency of meteorological data collected during 
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ABLE 2B can be found in Part I of this study and in both Garstang 
et al. (1990) and Harriss et al. (1990). To avoid repetition, 
only a concise summary of the data collection and measurement 

networks is presented below. 

The ABLE 2B measurement network was designed to capture 
atmospheric motions ranging from the planetary through the 
turbulent scale. Geostationary Observational Environmental 
Satellite (GOES) imagery was used to analyze cloud fields over 
the entire Amazon Basin and western tropical Atlantic Ocean. 

Hourly visible and infrared images with both a 2 km and 8 km grid 
resolution were also employed in partitioning the ACSL into 
convective end stratiform components. 

A four station mesoscale triangle network located in the 
CAB, 17 km north of the city of Manaus, Brazil (Fig. 1) was 
comprised of rawinsonde stations at the 3 triangle vertices co- 
located with Portable Automated Mesonet (PAM) stations. A PAM 
station was also employed at the center of the mesoscale 
triangle. Rawinsondes were launched every 3 h and provided wind, 
temperature, and humidity measurements between the surface and 
100 mb. The PAM stations, located on towers 5 m above the forest 
canopy, recorded measurements of surface parameters such as 
rainfall, u, v, q, T and 0 C at a frequency of 1 Hz and stored 
them as 1 min averaged values. 

Another component of the ABLE 2B measurement network was a 3 
cm radar located at the Eduardo Gomez Airport between Manaus and 

Constant Altitude Planned Position 


the mesoscale network. 
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Indicator (CAPPI) scans of the entire network area were made 
every 5 min and were used in tracking the convective cells 
embedded within the ACSL and propagating through the mesoscale 

triangle . 

As presented in Part I, the radar, PAM, rawinsonde and 
satellite data were used collectively to monitor the ACSL, to 
discern the stage of life cycle as it passed through the network, 
and to partition the ACSL into leading edge convection and 
trailing stratiform components. Heat and moisture transports for 
the different life cycle stages and the various ACSL components 
will be computed in the following sections. 

3. Budget Equations 

Equations presented by Yanai et al. (1973) and used in 
studies by Thompson et al. (1979), Johnson and Young (1983) and 
Johnson (1984), among others, are applied in heat and moisture 
transport calculations for the mesoscale triangle. The equations 
for sensible heat source (Q, ) and latent heat sink (Q 2 ) following 
Yanai et al. (1973) and Johnson and Young (1983), respectively, 

are : 


Q t = 


and 


d S 

3 t 


Q, — L 


+ V • Vs + u) 


3t 


3s 

3p 


+ V-Vq + 


3q 

3p 


(1) 

( 2 ) 


where the symbols are defined as follows: 
Q } = apparent sensible heat source 
Qv = apparent moisture sink 
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= dry static energy (= C p T + gz) 

= specific heat at constant pressure 
= temperature 
= acceleration of gravity 
= geopotential height 
= specific humidity 
= horizontal wind velocity 

= vertical velocity in pressures coordinates (= dp/dt) 
= pressure 
= time 

L = latent heat of condensation 

= area average of a quantity 

From Yanai et al. (1973) and Thompson et al. (1979), the equation 
tor the vertical eddy flux of total heat (F) between the 
tropopause <p,) and a given pressure level (p) is: 


s 

C P 

T 

g 

z 

q 

v 

til 

p 

t 


= g ' 1 


(Q, - Q 2 - Qr) d P 


(3) 


P, 


where Q„ is the net radiative heating rate and F is zero at the 
highest level. The right side of Eq. (3) is a good estimate for 
F in larger-scale studies such as Yanai et al. (1973) where the 
horizontal eddy flux convergences are considered negligible. 
However, as noted by Gallus and Johnson (1991), subgrid 
horizontal heat and moisture transports between the 

and stratiform components of a squall line may not be negligible 
on the mesoscale. 
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Thus, our estimate of F is actually an estimate of the total 
cloud transport, including both horizontal and vertical eddy 
fluxes. Although the vertical eddy flux term tends to dominate 
on the scale and size of systems we are dealing with, our values 
of F are still an overestimate of the vertical eddy flux of total 
heat and will be treated as such in later sections. 

Calculations of Q, and Q 2 include storage (3/3t), horizontal 
advection (V-V) and vertical advection (oj 3/ 3p) of s and q, 
respectively. Storage calculations use forward-in-time 
differences while vertical advection calculations use centered- 
in-space differences. Horizontal advection is found by applying 
a least squares objective analysis to rawinsonde measurements of 
s, q, u and v for the network triangle (Thompson et al., 1979; 
Johnson and Young, 1983) . 

The relatively small size of the triangle network and well- 
known limitations of rawinsonde systems cause unavoidable errors 
in Qi and Q,. These errors are smaller and less significant 
during convectively active periods (Thompson et al., 1979; 
Johnson and Young, 1983) when the vertical advection terms 
increase by an order of magnitude and dominate over storage and 

horizontal advection. 

Tropical budget studies (e.g., Thompson et al., 1979; Chen, 
1987) and methods of convective parameterization based upon them 
(Frank, 1978; Betts, 1986) have emphasized diabatic processes 
related to water, water vapor, and ice phase changes over 
radiation. Although the magnitude of Q R is small in relation to 
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Q, and Q„ Q r can play an important role in cloud development aa 
shown by model simulations (Tao et al., 1992), and can account 
for a 10-20% increase in surface precipitation. Studies by 
Webster and Stephens (1980), Gray and Jacobson (1977) and Tao et 
al. (1991) indicate that large differences in net radiation 
between the base and top of anvil clouds cannot be ignored in the 
total heat budget of tropical cloud systems. In addition, 
horizontal gradients resulting from radiation effects between 
clear and cloudy regions are also important. 

The net radiative heating (Qr) due to both longwave and 
shortwave radiative processes in the triangular volume is 
calculated by the methods described in detail by Pielke ( 

Separate equations are used for the calculation of radiative 
fluxes, both long and short wave, in clear air and m the 
presence of clouds . 

Fractional cloud cover (a,) at each pressure level in the 
mesoscale network is determined from the area-averaged relative 
humidity measured by rawinsonde. The method is generally based 
on cloud analysis procedures discussed in Air Weather Service 
Manual 105-124 (1969). The cloud analysis procedure employs a 

linear relationship between cloud cover and relative humidity. 

It is assumed that clouds form at a relative humidity of 90% and 
that the triangle is totally covered by clouds at a relative 
humidity of 100%. Relative humidities above the freezing level 
are based upon the frost point. The cloud analysis is checked 
for consistency against surface observations 


and satellite 
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images . 

For levels partially covered by clouds, the longwave 
radiative flux for the combination of clouds and cloud-free air 

is 

r 1w = a c R c + (1 - tf c )R (4) 

while the shortwave irradiance can be represented by 

r swt = a. R swc + (1 “ o : )K sw • (5) 

These methods neglect backscattered radiation from the ground and 
overlying clouds, cloud shading and the effects of cloud shape. 

Divergence for the mesoscale triangle is computed from the 
line integral of the wind components normal to each leg (Gamache 
and Houze, 1982) . The divergence is calculated from 1000 to 100 
mb in 25 mb intervals. The vertical velocity at each pressure 
level is first determined from the vertically integrated 
divergence and then adjusted by O'Brien s (1970) method. 

4. Results 

a. Heat and Moisture Transports 

Recent work on the total heat budgets of squall lines has 
shown the importance of partitioning the vertical heating and 
moisture distributions into convective and stratiform components 
(Johnson and Young, 1983; Johnson, 1984; Houze, 1989, Chong and 
Hauser, 1990; Gallus and Johnson, 1991). The rawinsonde network 
of the ABLE 2B experiment is smaller than either the convective 
or stratiform components of the ACSL. Thus, the network is 
capable of giving measurements (surface and upper air) obtained 
entirely within each component. This enables individual 
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calculations of heating and moistening profiles for both the 
convective and stratiform components. Throughout this study, a 
composite sounding derived from an average of the nearly 
simultaneously launched rawinsondes from the three corners of the 
mesoscale triangle (Fig. 1) will be used to compute heating and 
moisture profiles. 

Profiles of total sensible heat source and total apparent 
moisture sink for the 1 May 1987 squall line are shown in Fig. 2. 
These budgets were calculated for the system while it was in its 
mature stage. It should be noted that the profiles and budgets 
presented here for the 1 May system correspond only to that 
portion of the squall line that traversed the measurement network 
during the mature stage and may not be representative of other 

regions during different life cycles. 

The squall system produced net heating and drying throughout 
the troposphere with peaks at 550 and 500 mb, respectively. The 
shape of the heating profile closely resembles that diagnosed for 
mesoscale systems over the Gulf of Carpentaria, which also 
exhibit their heating peak at 550 mb during maturity (Frank and 
McBride , 1989 ) . 

Unlike the total heating profiles of continental squall 
lines studied over West Africa (Chong and Hauser, 1990) and the 
Great Plains (Callus and Johnson, 1991), the 1 May total Q, 
profile shows no cooling in the low-levels. Such cooling has 
been attributed to pronounced evaporation of both stratiform and 
convective rainwater in a relatively deep, dry subcloud air 
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layer. However, the lowermost troposphere over the Amazon forest 
canopy is characterized by a strong diurnal cycle of 
evapotranspiration that produces very large mixing ratios. 
Furthermore, model results for a similar ACSL system on 6 May 
( Scala et al., 1991) suggest a relative absence of dry 
tropospheric air for the production and maintenance of 
evaporatively cooled downdrafts. Thus, the characteristically 
high moisture content in the wet season lower troposphere may 
prevent or severely limit cooling in the convective region of the 

Amazon systems. 

Normalizing the convective apparent sensible heat source by 
rainfall rate allows comparison of the magnitude of the peak 
heating rate with other studies (Johnson, 1984) . The rainfall 
rate used in the normalization is expressed in cm d' 1 . The 
normalization rate is obtained by extrapolating the convective 
rainfall rates in mm min 1 to 24 h. A rainfall rate of 16 cm d' 1 
on 1 May gives a normalized value for the Ch maximum of 17.5°K d 
at 1200 UTC . This value is much larger than earlier estimates 
(Chong and Hauser, 1990; Gallus and Johnson, 1991; Johnson, 1976; 
Gamache and Houze, 1985). However, this large value may be 
partially due to the small size of the sampling network (1000 
km 2 ) relative to other studies. 

The profile of total apparent moisture sink in Fig. 2 
indicates a weak double peak structure characteristic of many 
tropical mesoscale convective systems (Johnson, 1984; Frank and 
McBride, 1989; Gallus and Johnson, 1991). In these systems, the 
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double peak in Q 2 results from the combined effects of a low- 
level peak in drying in the convective region and a midlevel 
drying maximum in the stratiform region (Johnson, 1984). 

However, Fig. 3, which partitions the total heating and drying 
profiles into convective and stratiform components, indicates 
that the peak drying is concentrated in the midlevels (450-650 
mb) in both the convective and stratiform regions. The weak 
double peak structure observed in the total profile is actually 
due to twin peaks (one at 500 mb, the other at 650 mb) in the 
convective region. This structure was also found m the deep 
convective elements over Taiwan (Johnson and Bresch, 1991). As 

noted in Johnson and Bresh (1991) , there are two possible 
explanations. The double peak may arise from either the strong 
convergence of eddy moisture flux in the deep convective towers 
and/or the mesoscale triangle sampled co-existing convective and 
stratiform cloud types within the budget volume. We are unable 
to distinguish between these two processes from the available 
data . 

Figure 3 also shows that the centers of convective heating 
and drying occupy nearly the same level and are approximately 
equal in magnitude. This suggests a coupling of Q, and Q 2 in the 
absence of pronounced eddy transports and is an indicator of 
predominantly stratiform processes (Gallus and Johnson, 1991) . 
study of systems over the Gulf of Carpentaria (Frank and McBride, 
1989) showed that, during the period of most intense rainfall 
(mature phase), total heating and drying were also coincident. 
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The stratiform region profile, Q 2 , in Fig. 3 exhibits two 
drying peaks, one at 800 mb and the other at 500 mb. While 
heating coincides with the upper peak as expected, cooling is 
associated with the lower maximum. In Part I, subsidence was 
found in this cooling layer. Typically, cooling and moistening 
characterize the sinking region beneath the melting level of most 
MCS (Houze, 1989). The reason for the observed drying is 
uncertain, but may be tied to the cutoff of a moisture source 
leading to drying of mesoscale downdrafts. 

The stratiform region also has a shallow layer of warming 
and moistening below 950 mb. Since the squall line passed 
through the sampling network very early in the morning, the 
warming and moistening may indicate rapid boundary layer recovery 
processes during the mid-afternoon solar heating of the forest 
canopy. These recoveries are indeed present, as indicated by the 
PAM surface time series of such values as q and 0 C (see Fig. 6, 
Part I) . Nearly instantaneous increases in sensible and latent 
heat fluxes are often measured just above the forest canopy in 
the wake of the squall systems. Formation of extensive stratus 
decks within and above the forest canopy immediately after the 
convective rain suggests that very large rates of 
evapotr ansp irat. ion are occurring. 

Figure 4 shows the vertical profile of the total cloud 
transport flux of heat (F) and its decomposition into convective 
and stratiform components. Once again, we acknowledge that our 
results for F are indicative of total cloud transport (both 
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horizontal and vertical) and are not indicative of the vertical 
eddy flux of total heat as calculated in larger scale studies 
(i.e., Yanai et al., 1973). However, in the convection region 
(LEC) , our F may still be a good estimate of vertical transports. 
In the stratiform region (TSR) though, horizontal motions are 
much more important and F overestimates the vertical eddy flux. 

The total profile shows a broad raid-tropospheric maximum 
which is composed of two distinct peaks, one at 450 mb and the 
other at 600 mb. These peaks are the result of the strong maxima 
in the convective eddy flux. The reason for the double 
convective maxima is unclear. The stratiform contribution to the 
total eddy flux is largest in the mid-troposphere, and is 
maximized at the same level as the convective flux. The large 
positive eddy surface layer value is associated with fluxes of 
moisture and heat during boundary layer recovery. Slightly 
negative values above this surface layer in the stratiform region 
may be a consequence of the downward transport of warm, dry air 
from an inversion layer between 700 and 950 mb (refer to Part I 
for the sounding structure) . 

The magnitude of the maximum total cloud transport at 600 mb 
is over 2200 W m' : , a number that nearly equals the maximum 
vertical eddy heat flux diagnosed between 500 and 650 mb for mia- 
latitude continental systems (Gallus and Johnson, 1991) . 

However, the mid-latitude total eddy flux does not exhibit a 
double peak structure. The values reported here are much larger 
than the 400 W m' 2 midlevel values found for tropical maritime 
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systems (Yanai et al., 1973; Esbensen et al., 1988), but it is 
important to note that in these earlier calculations, these 
budgets were calculated over much larger space and time scales 
than used in this study. 

Figure 5 presents the total cloud transport for the 
stratiform regions of three Amazon squall lines. The 1 May case 
represents a system in its maturity, the 26 April system was 
weakening and the 6 May case was a system that initially weakened 
but then began to regenerate. The maxima for all three systems 
has a similar magnitude, but each differs from the other in its 
vertical distribution. The 6 May regenerating case includes a 
secondary maximum located just above the surface at 950 mb. This 
may be associated with shallow (newly developing) convective 
processes as the system enters its reorganizing stage. 

Although the heat transport in the convective region of the 
1 May ACSL is larger and occupies a deeper layer than in the 
stratiform region, it will be shown that when integrated over the 
larger horizontal area of anvil cloud, the stratiform region 
makes a significant contribution to the vertical eddy heat flux 
of the total squall system, 
b . Svstem-wide transports 

The total system-wide heat transport of the ACSL is 
evaluated as a function of cloud component and life cycle with a 
model developed from the triangle network heat transport 
calculations and satellite analysis of active cloud area. 

The term "component" refers to part of the ACSL having features 
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that can be resolved by ABLE 2b sampling and analysis methods. 

As shown in Part I of this study, both the cloud component and 
life cycle stages of the ACSL as it passed through the rawinsonde 
measurement network were discerned through analysis of both GOES- 
west satellite imagery and the surface PAM data. 

We will show that vertical heat transport by ACSL can equal 
a substantial amount of the lateral heat export requirement for 
the tropics as estimated by Riehl and Simpson (1979). 

Furthermore, it is revealed that the dynamically active ACSL 
anvil plays an important role in the total system-wide heat 
transport . 

Table 1 gives the total vertical heat transport per unit 
area in both the LEC and TSR components of the mature 1 May ACSL. 
These values were obtained by integrating F over three separate 
layers. The integrations were carried out for the entire 
tropospheric column (1000-100 mb) , the upper troposphere (500 100 

mb) , and the layer of strongest Hadley export from the equatorial 
trough zone (300-100 mb) (Palmen and Vuorela, 1963) . 

As mentioned in Section 3, our estimate of F includes both 
the vertical and horizontal eddy fluxes of heat and is thus an 
overestimate of the vertical transport of total heat. The 
improved estimates of F, called F, and F 2 , have been made. In the 
Fj approximation the horizontal flux term is considered 
negligible. In the F ? approximation the horizontal flux is set 
equal in magnitude to the vertical eddy flux. Vertical motions 
dominate in the LEC, and therefore the true F will lie between F, 
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and F-, . In the TSR , where horizontal motions may dominate (see 
p ar t i) : the horizontal eddy flux term is (1) set equal to the 

vertical eddy flux (F t ) and (2) twice the vertical eddy flux term 
(Ft) . The vertical eddy flux of total heat (F) for both the 
convective (LEC) and stratiform (TSR) likely lies between F, and 

Ft . 

Although the total vertical heat transport within the LEC is 
2-3 times larger than in the TSR region, massive amounts of heat 
are transported upwards through the entire tropospheric column in 
both components. There is a substantial decrease in the vertical 
transport of heat in the 500-100 mb layer within the LEC. This 
decrease shows that the ACSL converts much of the energy 
originating in the lower levels to help maintain the convective 
towers and also to export hydrometeors from the LEC into the TSR 
where the heat is released. The heat reaching the 300-100 mb 
layer is only 2% and 6% of the column total in the LEC and TSR, 
but, as will be shown below, takes on global importance when 
integrated over the area of an entire ACSL. 

Table 2 displays the vertical heat transport per unit area 
in the TSR of a mature (1 May), weakening (26 April) and 
regenerating (6 May) ACSL. The general result of Tables 1 and 2 
for the ACSL anvil is unmistakable: it is dynamically active. 

Riehl and Simpson (1979) concluded that tropical deep 
convective towers provide the vertical heat transport required 
for the heat balance of the equatorial trough zone. In their 
previous study, Riehl and Malkus (1953), argued for the existence 
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of a midlevel total heat minimum, found throughout the tropics, 
that prohibited the vertical heat transport via a steady, large- 
scale upward circulation like that envisioned by Hadley (1735). 

Riehl and Malkus (1958) also concluded that undilute 
transport in the "hot tower" conduits of tropical thunderstorms 
is the necessary mechanism for the transfer of heat to the upper 
levels of the atmosphere. Budget studies by Johnson (1984) and 
model experiments by Scala et al. (1990, 1991, 1992) show that 
undilute tranpsort is not the only means of vertical heat 
transport in tropical squall lines over continental rainforests. 

Gamache and Houze (1982) and Cohen and Frank (1987) have 
shown that system-wide vertical transports are shared by cloud 
components consisting of more than hot towers. The results in 
Table 2 substantiate Johnson's (1984) recommendation to partition 
heating and moistening between the hot tower and anvil 
components . 

Substantial vertical heat transports occur in the TSR during 
all stages of the ACSL life cycle (Table 2) . The largest total 
column (1000-100 mb) heat transports occurred during the 
regenerating stage of the 6 May ACSL. However, the upper level 
vertical transports of heat (500-100, 300-100 mb) are strongest 

in the mature 1 May ACSL. This is especially true in the upper 
troposphere (300-100 mb) and is a good indicator of the 
importance of vertical transports in the stratiform region of a 
mature tropical squall line found over the CAB. 

The results of Tables 1 and 2 are extrapolated to the scale 
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of an entire ACSL with the aid of the Mb-enhanced satellite 
images. First, the length and width of each convectively active 
cloud cluster appearing in black (segment 7) is measured on the 
Mb images. Then the lengths and widths of all active clusters 
are summed over the length of the entire ACSL to obtain the total 

active cloud area for the system. 

The width of the LEC was derived from ground observations 

and the measured duration of convective rainfall at the PAM 
stations. The convective-stratiform rain rate cutoff was set at 
.25 mm/min. This value is slightly higher than most estimates 
(see for example, Gamache and Houze, 1982), but was necessitated 
by the .25 mm resolution of the tipping bucket raingauges 
employed at the PAM stations. However, the error introduced by 
this over-estimated cutoff is probably negligible due to the 
abrupt drop-off in rainfall intensities accompanying the passage 
of an ACSL (see Fig. 6, Part I) and the fast response rainfall 
measurements utilized on the PAM stations (every minute) . The 
widths obtained by these procedures are consistent with the 
Gamache and Houze (1982) structural morphology for tropical 
squall lines. 

The anvil width is obtained by subtracting the combined 
widths of the building cumulus and hot towers from the satellite 
measurement of the total active width of the system. Large 
differences in anvil width at different stages of system maturity 
result in large differences in total active cloud area during the 
COS life cycle. In the three ACSL evaluated, each had an active 
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length of about 750 km while passing over the mesoscale network. 

Table 3 shows the horizontally integrated system-wide 
heating in the 26 April, 1 May and 6 May ACSL for the same layers 
appearing in Tables 1-2. In these tables, it is assumed that 
vertical heat transport per unit area is uniform throughout each 
cloud component. 

In the case of the mature 1 May ACSL, the importance of the 
terms of total and vertical heat transport is clearly 
evident. Due to both the larger width of the anvil (when 
compared to the LEC) and the existence of vertical motions 
(incorporated into the calculations of Q t , Q 2 , and F) within it, 

there is almost as much total vertical heat transport within the 
TSR as in the LEC region. In the upper half of the atmosphere 
(500-100 mb), however, the vertical transports are larger in the 
TSR than in the LEC. However, a large portion of the energy 
present in the TSR and available for vertical transport was 
probably generated in the LEC and transported into the TSR. This 
is consistent with results from other studies on the heat budgets 

of tropical squall lines (Tao et al., 1992). 

All three ACSL cases show remarkably similar values of total 
vertical heat transport attributed to the TSR. The results in 
Table 2 most probably reflect the fact that, although the size 
(width) of the trailing anvil may change dramatically during the 
ACSL life cycle evolution, the strength of the circulations and 
transports within the core of the TSR remain relatively 
unchanged. Thus, the vertical motions and heat transports within 
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the TSR of all 3 storms were representative of an ACSL with 
active stratiform layer(s). In addition, the similar heat 
transport values support the classification of each storm as an 
ACSL and underscore some of their common attributes. 

The maximum poleward-directed mass flow in the upper level 
leg of the Hadley circulation occurs near 200 mb (Palmen and 
Vuorela , 1963). Riehl and Simpson (1979) calculated an upper- 

level requirement of 6 . 6 x 10 15 W for the lateral heat export from 
a io°-wide belt adjoining the equatorial trough line. 

Table 4 presents the relative contributions to the Riehl- 
Simpson export requirements of the 300-100 mb vertical heat 
transports within the 1 May LEC and TSR components, and within 
the TSR component of both the 26 April and 6 May cases. The 
relative contributions of the pre-storm cumulus heat transports 
in the 6 May and 26 April ACSL were found to be negligible (< 1%) 
and thus the total 300-100 mb heat transport is seen as a sum of 

the contributions made by the TSR and LEC. 

In the case of the 1 May ACSL, the 300-100 mb vertical 
transport of a mature system totals between 10% and 20% of the 
tropical export requirement. The percentage of the total 
tropical heat export requirement accomplished by a regenerating 
(6 May) and weakening (26 April) ACSL is found to be 6-10% and 
10-16%, respectively. 

According to analyses of ABLE 2b satellite data, up to 3 
ACSL at different stages of their life cycle may exist 
concurrently over tropical South America. Using the results 
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obtained in Table 4, it can be seen that, during active periods 
of ACSL formation in the Amazon Basin (April-May and August- 
September) (Cohen et al., 1989), the total instantaneous upper- 
level heat transport attributed to the ACSL may equal as much as 
or higher than one-third of the tropical heat export requirement. 

Similar, although slightly smaller values were determined by 
Connors (personal communication) . In a study which classified 
the weather regimes over the Amazon Basin through analysis of 
digital data from hourly GOES-6 imagery (1 April 1986-1 April 
1987), Connors (personal communication) found several prominent 
wet season regimes. The wet season weather regimes which existed 
during the climatological ACSL active months of April-May (Cohen 
et al., 1989) often contributed 1. 0-2.0 x 10 15 W or 15-30% of the 
energy required to balance the heat budget of the global tropics. 

Riehl and Malkus (1958) argued that 1500-5000 hot towers, 
distributed among approximately 30 tropical synoptic disturbances 
could meet the heat export requirements of the equatorial trough 
at any one time. Our results indicate that 5-10 organized meso- 
to synoptic-scale disturbances such as the ACSL can meet the 
instantaneous upper-level heat transport requirement for the 
tropics . 

5. Summary and Conclusions 

In Part II of our investigation of the Amazon Coastal Squall 
Lines (ACSL) we have analyzed the heat and moisture transports of 
the ACSL in terms of convective and stratiform components, 
instantaneous and storm-wide energetics, and in relation to the 
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energy required to balance the heat budget of the global tropics. 

Using rawinsonde data collected from a 4-station mesoscale 
network (approximately 1000 km 2 ) located in the CAB during the 1 
April-15 May 1987 ABLE 2b experiment, vertical profiles of Q, , Q? 
and F were determined for the convective (LEC - Leading Edge 
Convection) and stratiform (TSR - Trailing Stratiform Region) 
regions of a mature ACSL (1 May 1987) and the stratiform region 
of a -weakening (26 April 1987) and regenerating (6 May 1987) 

ACSL . These instantaneous profiles were extrapolated to storm- 
wide energetics using rainfall rate and satellite imagery data to 
delineate the individual components of the ACSL and estimate 

their size and active cloud area. 

Analysis of the mature 1 May ACSL indicates net heating and 

drying through the troposphere. The peak heating was found in 
the mid-levels (550-500 mb). This is similar to results for 
mesoscale systems in the Gulf of Carpenteria (Frank and McBride, 
1989) but somewhat above the level of peak heating for mature 
eastern Atlantic systems (Johnson, 1984) and below the level 
identified for West African (Chong and Hauser, 1990), western 
Pacific (Johnson, 1984) and U.S. Great Plains (Gallus and 
Johnson , 1991 ) . 

Unlike the double peak structure found in the vertical 
profiles of Q 2 for many tropical mesoscale systems, the peak 
drying in the 1 May ACSL is mainly concentrated in mid-levels 
(450-650 mb) . The similarity in magnitude and level of maximum 
transport found in the profiles of total Q[ and Q 2 suggests a 
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coupling of Q, and Q 2 in the absence of pronounced eddy 
transports. This is also indicative of a predominantly 
stratiform process (Gallus and Johnson, 1991) . 

The Qj , Q, and F profiles for the TSR of the 1 May, 26 April 
and 6 May ACSL all displayed similar shapes and magnitudes, with 
peak values found between 550 and 700 mb. Although the 
instantaneous heat and moisture fluxes associated with the LEC 
were much stronger than in the TSR, the horizontal and vertical 
extent of the TSR over the entire storm system results in 

significant heat transports for the TSR. 

Total cloud transports of heat (F) were obtained for the 
atmospheric column (1000-100 mb) and the upper troposphere (500- 
100 mb, 300-100 mb) for both the LEC and TSR of the 1 May ACSL. 
Previous studies on larger scale networks (i.e., Yanai et al., 
1973) have considered horizontal transports negligible and thus F 
was a representation of total vertical eddy flux of heat. 

However, in the present mesoscale study, horizontal effects can 
not be neglected and thus F is a measure of the total cloud 
transport (horizontal and vertical) . Despite this difference, 
the similarity between our results and other studies indicate the 
critical role of cloud scale transports in the overall 
distributions of Q,, Q 2 and F. 

Large amounts of heat were transported in both the LEC and 
TSR components; 250,000-520,000 W m' in the LEC and 170,000 
340,000 W m' 2 in the TSR. The heat reaching the 300-100 mb layer 
was only 3% (LEC) of the total column heating. Thus, much of the 
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low-level energy is used to both maintain the convective towers 
and to export hydrometeors into the anvils. 

The TSR of the regenerating (6 May) and weakening (26 April) 
ACSL also produced substantial heat transport of over 200,000 W 
nr 2 . Thus, a dynamically active anvil region is a characteristic 
of the ACSL. In all three TSR, only 2-6% of the total column 
heat transport reaches the 300-100 mb layer for export from the 
Amazon Basin. However, calculations show that these small 
amounts are significant when integrated over the active area of 
the entire ACSL. 

ACSL lengths and. widths were determined through analysis of 
surface and satellite data. Typical dimensions of an active ACSL 
were 750-1000 km long by 100-200 km wide. In all 3 cases, the 
TSR occupied over 65% of the ACSL width. 

On a storm-wide basis, the TSR region of the 1 May ACSL was 
responsible for more total vertical heat transport than the LEC 
region. This is especially true in the 500-100 mb layer. 

Similar heat transports were also found in the TSR of the 26 
April and 6 May ACSL, underscoring their importance to the 
dynamics and energetics of the ACSL as a whole. 

The relative contribution to the global heat budget was 
determined. The upper levels of a single, mature ACSL 
could account for 15-20% of the total amount of heat export 
required to balance the heat budget of the global tropics (Riehl 
and Simpson, 1979). 

Although our calculations of F represented total cloud scale 
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transports, an estimate of the total vertical eddy flux of heat 
was obtained by assuming in the convective region (LEC) that, (1) 
as in larger scale studies (Gallus and Johnson, 1991; Yanai et 
al., 1973), horizontal transports are negligible; and (2) the 
vertical and horizontal transports are equal in magnitude. In 
the stratiform region (TSR) , the horizontal flux was (1) set 
equal to and (2) 2X larger than the vertical flux. The true 
value of F likely lies somewhere between the 2 values. 

Analysis of ABLE 2b satellite data indicates the presence of 
three ACSL over the Amazon Basin on several occasions. 
Calculations show that, in these instances, the ACSL may be 
responsible for up to one-third of the required upper level heat 
export from the global tropics. Thus, it appears that, in 
addition to the Riehl and Malkus (1958) hot tower hypothesis, 5- 
10 organized meso- to synoptic-scale disturbances (such as the 
ACSL) can meet the instantaneous heat export requirement for the 


tropics . 
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TABLE 1. Vertical heat transport (per unit area) in the 
mature 1 May 1987 ACSL as a function of cloud component. 


Pressure Layer 
(mb) 


LEC 

(W m : ) 


TSR 
(W m' 2 ) 


1000-100 

500-100 


521, 538 
112 , 350 


260,769 
56 , 175 


171,588 
64 , 690 


114 , 380 
43 , 122 


300-100 


12 , 950 


6,475 


9 , 299 


6, 198 



TABLE 2. Vertical heat transports (per unit area) in 
the TSR region of a weakening (W) , mature (M) and 
regenerating (R) ACSL. 



1000- 

(W 

F, 

-100 mb 

m 2 ) 

F, 

500-100 mb 
(W ra' 2 ) 

F, F 3 

300-100 mb 
(W m' 2 ) 

F, 

26 April (W) 

182 , 700 

121,787 

51,671 

34 , 443 

8 , 297 

5,531 

1 May (M) 

171, 588 

114,380 

64 , 690 

43 , 122 

9 , 299 

6 , 198 

6 May (R) 

206,294 

137,516 

57 , 167 

38 , 108 

6 , 975 

4 , 650 



TABLE 3. Horizontal dimensions and total vertical heat transport 
in 3 separate layers for a mature (1 May), weakening (26 April) and 
regenerating (6 May) ACSL. Definitions for F, and F 2 are provided 
in text (Section 4b) . Heat transports (W) were computed for the 
cloud component sampled by the rawinsondes. 

Mature : 1 MAY 



PSC 

LEC 

F, F 2 

F, 

TSR 

F 2 

Width (km) 

25 

o 

in 


100 


Length (km) 

750 

750 


750 


1000-100 mb 

- 

2 . 0xl0 16 W 1.0X10 16 

1. 3xl0 16 


0 . 9xl0 16 

500-100 mb 

- 

4 . 2xl0 15 W 2 . lxlO 15 

4 . 8xl0 15 


3 . 2xl0 13 

300-100 mb 

- 

4 . 8xl0 14 W 2.4X10 14 

7 . OxlO 14 


4 . 6xl0 14 


Weak: 26 APRIL 



PSC 

F t F 2 

LEC 

F, 

TSR 

f 3 

Width (km) 

25 

50 


100 

Length (km) 

750 

i 750 


750 

1000-100 mb 

5 . 2xlO l5 W 2 . 6xl0 15 

- 

1. 8xl0 16 

0 . 9xl0 16 

500-100 mb 

3 . 6xl0 14 W 1 . 8x1 0 14 

- 

3 . 8xl0 15 

2 . 5xl0 13 

300-100 mb 

2 . 6x 10 13 W 1 . 3xl0 13 

- 

6 , 2xl0 14 

3 . 9xl0 14 



Table 3 . (cont . ) 



TABLE 4 . 


The 300-100 mb ACSL heat transport in percent of the 
equatorial heat export requirement of 6.6 x 10 15 W (Riehl and 
Simpson, 1979) . See Table 3 for definition of F, and F 2 . 



LEC 

(%) 

F, F 2 

TSR 

(%) 

Fi F 2 

TOTAL 

(%) 

F, F 2 

1 May (mature) 

7 . 3 

3 . 6 

10 . 6 

7 . 0 

17 . 9 

10 . 6 

26 April (weakening) 

*7 . 3 

3 . 6 

9 . 0 

5.9 

*16.3 

9 . 5 

6 May (regenerating) 

*3 . 7 

1 * 9 

5 . 9 

3 . 9 

* 9.6 

5 • 8 


* Assuming the LEC 300-100 mb transport per unit area to be no 
larger than that occurring in the mature 1 May LEC region. 
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